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Chapter  1 
Introduction 


Since  its  synthesis  in  the  1920s  by  Goldschmidt  [1],  Gallium  Arsenide  has  received 
much  attention  in  the  last  few  decades.  In  the  mid-1980s,  GaAs  technology  finally 
matured  into  the  age  of  production.  We  saw  a  boom  of  companies  dedicated  to  the 
growth  of  GaAs  materials  and  the  fabrication  of  GaAs  devices  and  integrated  circuits. 
Although  GaAs  is  no  longer  being  considered  a  general  purpose  material  like  silicon, 
it  is  now  well  established  in  several  niche  markets,  such  as  Direct  Broadcast  Satellite, 
Microwave  Monolithic  Integrated  Circuits  and  Optoelectronics. 


1.1  Physical  Characteristics  of  GaAs 

Gallium  arsenide  is  composed  of  the  group  III  element  Ga  and  the  group  V  element 
As.  Each  of  these  two  elements  form  a  sublattice  of  a  face-centered-cubic  structure. 
Figure  1.1  shows  the  Ga  and  As  arrangement  in  a  conventional  unit  cell.  The  Ga 
and  As  sublattices  are  offset  by  [1—j.  This  crystal  configuration  is  known  as  cubic 
sphalerite  or  zincblende.  Therefore,  the  bond  length  between  nearest  neighbors  of  Ga 
and  As  atoms  is  given  by,  r0  =  >/3A/4,  where  A  is  the  length  of  the  unit  cube. 

The  size  of  the  unit  cell  for  undoped  GaAs,  both  melt-grown  and  epitaxial  at 
300K,  is  given  by  [2]: 

A 300  =  5.65325  ±  0.00002A  (1.1) 
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However,  the  GaAs  lattice  constant  can  increase  by  up  to  0.02%  for  large  concen¬ 
trations  of  dopants,  such  as  Sn  or  Te.  In  addition,  non-stoichiometric  GaAs  has  a 

slightly  smaller  value  of  A  for  the  As-rich  crystal  and  slightly  greater  A  for  Ga-rich 
crystal. 

GaAs  also  has  the  natural  cleavage  plane  of  (110).  Figure  1.2  shows  a  unit  cell 
bisected  by  such  a  plane.  This  type  of  plane  contains  both  Ga  and  As  atoms.  Another 
important  plane  in  GaAs  is  the  (111)  plane,  shown  in  Figure  1.3.  This  family  of  planes 
only  contains  one  type  of  atom.  As  a  result,  (111)  and  (III)  planes  are  not  equivalent. 
Those  that  only  contain  Ga  atoms  are  often  called  {111 A)  and  those  with  only  As 
atoms  are  called  {111B}.  This  distinction  is  important  because  the  different  planes 
have  different  chemical  properties;  therefore,  have  different  growth  and  etching  rates. 
Table  1.1  summarizes  the  important  properties  of  GaAs  at  300K. 


s 


.2:  The  (110)  Plane  of  GaAs. 


.3:  The 


(111)  Plane  of  GaAs. 


CHAPTER  1.  INTRODUCTION 


Molecular  Weight  (M) 

144.6416 

Lattice  Constant,  A 

5.65325 

Density  of  Solid,  d,  (g  cm"3) 

5.307 

Density  of  Liquid,  d\  (g  cm-3) 

5.71 

Atomic  Concentration,  Na  (cm-3) 

4.42  x  1022 

Melting  Temperature,  Tm(C) 

1238  -♦  1240 

Vapor  Pressure  at  Tm,  pm  (bar) 

0.976 

Thermal  DifFusivity,  Dt  (cm2s-1) 

0.218 

Thermal  Conductivity,  K  (Wcm-1K-1) 

0.49 

Linear  Expansion  Coefficient,  a  (K-1) 

5.73  x  10-6 

Specific  Heat,  (Jg-1K-1) 

0.35 

Molar  Heat  Capacity,  Cp  (Jmol-1K_1) 

47.026 

Standard  Enthalpy,  Ho  (kJmol-1) 

-81.643 

Standard  Entropy,  S0  (Jmol-1K-1) 

64.267 

Standard  Free  Enthalpy,  Go  (kJmol-1) 

-100.806 

Enthalpy  for  Melting,  LHm  (kJmol-1) 

96.7 

Stacking  Fault  Energy,  E,f  (mJm-2) 

55.0 

Hardness  (on  Mohs  scale) 

4  — >  5 

Surface  Microhardness  (Knoops  Test  kgmm-2) 

750  ±  40 

Table  1.1:  Physical  Properties  of  GaAs  at  300K 
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1.2  Crystal  Growth  and  Processing  Techniques 

GaAs  is  grown  mainly  by  two  methods,  horizontal  Bridgman  (HB)  and  Liquid  en¬ 
capsulated  Czochralski  (LEC).  Most  commercial  wafers  today  are  grown  by  the  later 
method,  although  HB  grown  wafers  are  still  being  used  in  research. 

1.2.1  Horizontal  Bridgman  (HB) 

The  starting  material  for  the  HB  method  is  pure  gallium  or  polycrystalline  GaAs, 
placed  in  a  long  boat.  The  boat  is  sealed  in  a  quartz  ampoule  which  is  typically  filled 
with  an  inert  gas.  A  GaAs  seed  crystal  is  placed  in  the  front  of  the  boat.  Elemental 
arsenic  is  also  present  in  the  ampoule.  This  configuration  is  shown  in  Figure  1.4. 

The  ampoule  is  placed  in  a  furnace  made  up  of  a  number  of  consecutive  heating 
zones.  The  temperature  along  the  length  of  the  ampoule  is  carefully  regulated  to 
achieve  the  desired  growth.  The  elemental  arsenic  is  brought  to  the  gaseous  state, 
while  the  Ga  or  poly  GaAs  is  heated  to  the  melting  point  of  GaAs.  The  temperature 
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of  the  seed  region  is  just  below  the  GaAs  solidification  point  [1].  The  crystal  is  usually 
grown  in  the  <  111  >  direction,  which  is  the  slow  growth  direction. 

The  completed  crystal  has  a  shape  that  matches  that  of  the  boat.  When  (100) 
oriented  wafers  are  desired,  HB  crystals  are  cut  at  an  angle  of  54.7°  to  the  ingot  axis. 
By  sawing  the  ingot  at  an  angle,  the  compositional  variations  along  the  axis  of  the 
crystal  becomes  a  variation  across  the  wafer;  therefore,  HB  wafers  are  undesirable  for 
IC  type  applications. 

1.2.2  Liquid  Encapsulated  Czochralski  (LEC) 

The  starting  point  for  the  LEC  method  is  molten  GaAs.  The  melt  is  confined  by 
a  layer  of  liquid  boric  oxide  (B2O3).  Both  the  melt  find  B2O3  are  contained  in  a 
crucible.  Figure  1.5  shows  the  setup  of  the  LEC  method.  Crystal  growth  is  started 
when  the  GaAs  seed  penetrates  the  boric  oxide  and  makes  contact  with  the  melt. 
The  diameter  of  the  crystal  is  controlled  by  the  pull  rate  [1]. 

The  crucible  used  in  the  LEC  method  is  either  quartz  (Si02)  or  pyrolytic  boron 
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nitride  (PBN).  The  quartz  crucible  adds  Si,  which  is  a  n-type  dopant,  to  the  grown 
crystal.  In  order  to  maintain  the  semi -insulating  property  of  the  GaAs  wafer,  Cr,  a 
deep  acceptor,  is  intentionally  added.  However,  Cr  doped  material  has  many  unde¬ 
sirable  effects  such  as  outdiffusion  during  heating  and  low  distribution  cofficient.  The 
use  of  quartz  crucibles  is  no  longer  widespread. 

GaAs  crystal  grown  with  pyrolytic  boron  crucibles  still  contain  carbon,  which 
comes  from  the  graphite  heaters.  Carbon  is  a  shallow  acceptor  and  it  is  compensated 
by  the  deep  donor  EL2  so  the  material  retains  the  desired  semi-insulating  property. 

1.3  Epitaxy 

Epitaxy  is  the  process  of  growing  GaAs  material  on  the  surface  of  a  GaAs  substrate. 
Ga  and  As  atoms  are  brought  to  the  GaAs  surface  under  suitable  temperature  and 
concentration.  The  epitaxial  layer  of  GaAs  is  typically  of  much  higher  quality  than 
the  substrate. 

In  molecular  beam  epitaxy,  the  substrate  is  placed  in  a  vacuum  and  elemental 
species  are  evaporated  from  ovens  and  on  to  the  heated  substrate.  Figure  1.6  illus¬ 
trates  the  schematics  of  a  molecular  beam  epitaxy  system.  The  advantage  of  MBE 
is  that  it  can  produce  almost  any  epitaxial  layer  composition,  thickness  and  dop¬ 
ing.  However,  the  system  requires  high  vacuum  setup  as  well  as  other  complex  and 
expensive  equipment.  The  growth  is  also  very  slow,  typically  around  1  /xm/h. 


1.4  Ion  Implantation 

Ion  implantation  at  present  is  the  most  economical  and  widely  used  method  for  pro¬ 
ducing  active  layers  for  GaAs  FET  or  MMIC  devices.  In  this  process,  ions  are  im¬ 
planted  into  the  substrate  at  an  energy  of  30  to  400  keV  at  doses  of  1012  to  1014 
atoms  per  cm2.  The  implantation  process  creates  damage  to  the  substrate  and  the 
wafers  must  be  subjected  to  a  high  temperature  annealing  step  to  repair  the  lattice 
damage  and  allow  the  implanted  atoms  to  move  into  the  proper  lattice  site.  This 
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Figure  1.6:  Molecular  Beam  Epitaxial  of  GaAs. 


step  is  known  as  activation  of  the  implant.  Not  all  dopants  however,  become  incorpo¬ 
rated  into  a  proper  lattice  site.  Therefore,  100%  activation  is  usually  achieved,  and 
is  typically  found  to  lie  in  the  range  of  70%  to  90%. 

The  major  advantages  and  disadvantages  of  ion  implantation  are  summarized  as 
follows:  [3] 

•  Advantage: 

1.  Reproducibility,  uniformity  and  high  speed  of  the  doping  process. 

2.  Exact  control  of  the  amount  of  dopant. 

3.  Less  stringent  requirements  on  dopant  source  purity  since  mass  separation 
is  used. 

4.  Avoidance  of  high  temperature  during  the  implant  process  itself. 

5.  Ability  to  use  simple  masking  methods  which  allows  selective  doping. 


•  Disadvantage: 
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1.  Creation  of  radiation  damage  which  requires  high  temperature  annealing 
for  repair. 

2.  Secondary  effects,  such  as  diffusion  during  annealing  or  channeling  during 
implantation  alter  the  expected  dopant  profile. 


Chapter  2 

Basic  Semiconductor  Device 
Physics 


The  set  of  partial  differential  equations  that  govern  potential  distribution,  carrier 
concentrations  and  current  flow  in  any  semiconductor  device  are  called  the  basic 
semiconductor  device  equations.  They  may  be  derived  from  the  Maxwell’s  equations 
2.1  to  2.4. 


dD 

(2.1) 

dB 

v*£  =  -ar 

(2.2) 

VD  =  p 

(2.3) 

V  •  £  =  0 

(2.4) 

where  E  and  D  are  the  electric  field  and  displacement  vectors  respectively.  H  and 
B  are  the  magnetic  field  and  induction  vectors.  The  conduction  current  density  is 
denoted  by  J  and  p  is  the  electric  charge  density. 
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2.1  Poisson’s  Equation 

To  derive  Poisson’s  equation,  we  need  to  first  introduce  a  relationship  between  D  and 
E.  This  is  given  by: 

D  =  eE  (2.5) 

where  £  is  the  permittivity  of  the  medium.  We  assume  here  the  permittivity  is 
time  independent  and  spatially  homogeneous  in  the  semiconductor.  In  general,  e  is 
represented  by  a  3  x  3  matrix.  We,  however,  assume  that  the  material  is  isotropic  and 
thus  e  is  given  by  a  scalar.  For  the  usual  applications  in  semiconductor  devices,  these 
assumptions  are  quite  realistic.  For  GaAs,  e  =  13.l£v  where  ev  is  the  permittivity 
constant  in  vacuum  [4]. 

Due  to  equation  2.4,  we  may  introduce  a  vector  potential  such  that: 

B  =  Vx  A,  V  •  A~0  (2.6) 
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finally  arrive  at  Poisson’s  equation: 


eV2V>  =  —p 


(2.11) 


For  semiconductor  device  applications,  the  space  charge  density  can  be  further 
broken  down  into  three  parts: 


p  =  q(p-n  +  C)  (2.12) 

where  p  is  the  concentration  of  positively  charged  holes  and  n  is  the  negatively  charged 
electron  density.  C  represents  the  ionized  dopant  profile  and  is  given  by  the  difference 
between  the  ionized  donor  and  acceptor  profiles. 


2.2  Continuity  Equations 


Applying  the  divergence  operator  to  the  first  Maxwell  equation  2.1  and  using  2.3: 

°  =  V-J  +  ^  (2.13) 


We  now  split  the  conduction  into  two  parts,  Jn  and  Jp,  resulting  from  electrons  and 
holes  respectively: 

J  =  Jn  +  J„  (2.14) 

We  also  assume  that  the  ionized  charge  profile  is  time  independent,  ie, 


(2.15) 


This  effectively  neglects  contributions  from  any  charged  defects,  such  as  vacancies, 
dislocations  and  deep  recombination  traps,  which  may  change  their  charge  state  over 
time. 

Then  using  2.14,  2.15  and  2.12,  we  obtain: 


T-r  t  t-r  r  dn 

-  v’Jp~qdi  ~  v'Jn~qm 


(2.16) 
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By  setting  the  left  and  right 


side  of  the  above  equation  equal  to  qR,  we  obtain: 


dn 

VJn~qdi  =  qR 


(2.17) 


V  •  Jp  +  =  -qR  (2.18) 

Obviously,  we  can  not  gain  any  additional  information  by  splitting  a  single  equation 
into  two  equations.  However,  these  formal  steps  can  help  us  interpret  the  equations 
more  easily.  The  quantity  R  may  be  understood  as  the  difference  of  the  rate  at  which 
electron  and  hole  carriers  pairs  recombine  and  the  rate  at  which  they  are  generated  in 
the  semiconductor.  Therefore,  we  call  R  the  recombination-generation  rate.  Positive 
R  would  mean  net  recombination  and  negative  R  means  net  generation.  However,  if 
R  is  modeled  with  the  knowledge  of  solid  state  and  statistical  physics,  each  of  these 
two  equations  has  a  meaning  of  its  own  and  should  be  considered  as  two  equations. 


2.3  Current  Equations 

We  assume  here  there  are  two  types  of  current  in  the  semiconductor  device.  The 
diffusion  of  the  electron  and  hole  ensembles  results  in  diffusion  current  J*'? I  and 

the  drift  of  electrons  and  holes  caused  by  the  electric  field  results  in  drift  current  J 
Jpr,7‘.  We  also  assume  that  the  total  current  is  determined  by  linearly  superimposing 
the  diffusion  and  drift  currents: 


Jn  =  4iff  +  J?,f\  Jv  =  J*fi  +  Jfift  (2.19) 

Electron  and  holes  diffuse  from  regions  of  high  respective  concentrations  to  regions 
of  lower  concentrations.  The  direction  of  such  diffusion  is  the  direction  of  steepest 
descent  of  particle  concentration.  Therefore,  the  diffusion  flux  densities  must  be 
proportional  to  the  gradients  of  the  particle  concentration.  As  a  result,  the  diffusion 
currents  are  given  by: 


Jf "  =  qD „Vn 


(2.20) 
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J}if/  =  -qDpVp  (2.21) 

The  signs  are  chosen  here  so  that  the  diffusion  coefficients  Dn  and  Dp  are  positive. 
The  drift  currents  are  given  by: 


s. 

II 

(2.22) 

= wt 

(2.23) 

where  v*  and  vp  are  the  average  drift  velocity  of  electrons  and  holes  respectively.  The 
drift  of  holes  has  the  same  orientation  as  the  electric  field  while  the  drift  of  electrons 
has  the  opposite  direction.  The  drift  velocities  are  related  to  the  electric  field  by: 

vdn  =  - UnE ,  vdp  =  ppE  (2.24) 

where  pn  and  pp  are  the  electron  and  hole  mobilities.  By  inserting  equations  2.22, 
2.23  and  2.24  into  2.19,  we  obtain  the  current  relations: 

Jn  =  g£>„Vn  +  qp„nE  (2.25) 

Jp  =  -qDpVp  +  qpppE  (2.26) 

Clearly,  the  above  equations  are  derived  from  a  phenomenological  approach  instead 
of  well  established  physical  principles.  Therefore,  the  current  equations  are  not  to 
be  accepted  in  the  same  way  as  Poisson’s  equation  and  the  continuity  equations. 
Their  validity  must  be  investigated  by  other  means.  A  more  rigorous  approach  to  the 
derivation  of  the  current  equations  leads  to  the  same  expressions  by  starting  with  the 
Boltzman  transport  equation  [7]. 

2.4  Summary 

The  basic  semiconductor  device  equations  are  given  as: 


Poisson's  Equation  e Vfy  =  q(n  — p—C ) 


(2.27) 


4 
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Electron  Continuity 


V  Jn 


—  qR 


Hole  Continuity 
Electron  Current 
Hole  Current 


VJP  +  q-~  =  -qR 
Jn  =  qHnnE  +  $Z)„Vn 


Jp  =  qUppE  -  qDpVp 


(2.28) 

(2.29) 

(2.30) 

(2.31) 


2.5  Basic  Operation  of  GaAs  MEFSET 

Unlike  Si,  GaAs  does  not  have  a  native  oxide,  so  devices  such  as  the  Metal  Oxide 
Silicon  Field  Effect  Transistor  (MOSFET)  can  not  be  fabricated.  The  most  dominant 
device  in  both  digital  and  analog  applications  is  the  Metal  Semiconductor  Field  Effect 
Transistor  (MESFET).  It  closely  resembles  a  Junction  Field  Effect  Transistor  (JFET) 
but  uses  a  metal  gate  to  control  the  current  transport  in  the  channel  by  modulating 
the  extent  of  the  depletion  region  formed  by  the  Schottky  barrier. 

2.5.1  Metal  Semiconductor  Junction 

When  a  metal  and  a  semiconductor  are  brought  into  intimate  contact,  electron  trans¬ 
fer  will  take  place  between  the  systems.  When  equilibrium  is  reached,  every  process 
and  its  reverse  is  occuring  at  the  same  rate  and  no  net  electron  transfer  occurs.  The 
two  system  will  thus  have  the  same  Fermi  energy. 

Let’s  define  the  work  function  as  the  energy  difference  between  the  Fermi  and  the 
vacuum  levels.  Let  q<f>m  be  the  work  function  of  the  metal  and  q(\  +  Vn)  be  the  work 
function  for  the  semiconductor.  Here  \  is  the  electron  affinity  measured  from  the 
conduction  band  Ec  to  the  vacuum  level  and  V„  is  the  energy  difference  between  Ec 
and  Fermi  level.  These  definitions  are  illustrated  in  Figure  2.1. 

As  the  contact  is  formed,  electrons  at  the  Fermi  level  in  the  metal  see  a  barrier 
°f  q<f>Bn  and  the  electrons  at  the  conduction  band  edge  in  the  semiconductor  see  a 
barrier  height  of  qVei.  In  equilibrium,  we  have  [6]: 


q<l>Bn  =  ~  X) 


(2.32) 
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VACUUM 


Figure  2.1:  Ideal  Metal  Semiconductor  Junction. 


qVBi  =  q(<t>m-x)-qVn  (2.33) 

The  above  equations  only  hold  if  there  are  no  surface  states.  Usually,  however,  semi¬ 
conductors  have  significant  amount  of  surface  states  which  often  lie  in  the  bandgap. 
Therefore,  real  metal  semiconductor  junctions  are  never  quite  in  accord  with  the  above 
model.  Si  comes  very  close  to  the  ideal  situation,  but  most  compound  semiconductors 
do  not  [7]. 

In  the  case  of  large  numbers  of  surface  states,  as  shown  in  Figure  2.2,  all  the 
electron  transfer  required  to  achieve  equilibrium  can  be  accomodated  at  the  surface. 
Therefore,  space  charge  in  the  semiconductor  remains  unaffected.  As  a  result,  the 
barrier  height  is  determined  by  the  properties  of  the  semiconductor  alone  and  is 
independent  of  the  metal.  This  case  applies  to  semi-insulating  GaAs  and  the  experi¬ 
mentally  measured  barrier  height  is  found  to  be  0.78  eV. 

The  charge  distribution  diagram  for  an  ideal  metal  semiconductor  junction  is 
shown  in  Figure  2.3.  Since  metal  closely  approximates  a  perfect  conductor,  the  charge 
transfered  to  it  from  the  semiconductor  exists  only  on  the  surface  in  a  thin  plane.  In 
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VACUUM 


Figure  2.2:  Metal  Semiconductor  Junction  with  Large  Surface  State. 


the  n-  type  semiconductor,  the  charge  consists  of  ionized  donors,  electrons  and  holes. 
The  hole  concentration  is  extremely  small  since  it’s  the  minority  carrier.  Under  the 
depletion  approximation,  we  also  assume  that  the  electron  concentration  in  the  deple¬ 
tion  region  which  extends  to  x=xj  is  much  less  than  the  ionized  donor  concentration. 
Outside  of  the  depletion  region,  the  electron  and  ionized  donor  densities  are  assumed 
to  be  equal. 

Under  the  depletion  approximation,  the  Poisson’s  equation  may  be  transformed 
as: 

-  f?  =  %  *  !  fa*)  -  »<*)  +  A®(x)  -  N;(z))  (2.34) 

d7V  q 

— xt  ~  -No  for  0  <  x  <  Xd  (2.35) 

ox2  e 

The  electric  field  €  can  then  be  solved  for  as: 

€{x)  =  Em-  =  S**D{Xi  -  x)  (2.36) 


where  Sm  is  the  maximum  field  that  exists  at  x=0.  The  potential  may  be  obtained 
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Figure  2.3:  Charge  Distribution  in  Metal  Semiconductor  Junction. 


by  integrating  the  electric  field  equation. 

V(x)  =  ^-(xix-\z2)  +  4’B  (2.37) 

The  depletion  width  may  be  obtained  by  setting  x  =  xj  in  the  above  equation  where 
the  potential  is  equal  to  Vn  and  H,  =  <f>s  —  Vn.  So  we  obtain: 

IJ  =  '  V))  (2-38) 
where  V  is  the  applied  potential. 

2.5.2  Constant  Mobility  Model  of  MESFET 

The  assumption  of  constant  mobility  is  accurate  in  the  case  of  long  gate-length  FET 
where  the  field  is  sufficiently  small.  Figure  2.4  illustrates  the  geometry  of  such  a 
device.  The  source  here  is  grounded  and  drain  biased  to  a  small  positive  voltage. 


CHAPTER  2.  BASIC  SEMICONDUCTOR  DEVICE  PHYSICS 


19 


z=0  z=L 

< - L - ► 


Source  Gate  Drain 


The  gate  is  negatively  biased  and  the  absolute  value  of  the  gate  bias  is  VG.  Then,  the 
depletion  layer  width  h(z)  varies  slowly  along  the  channel.  The  electric  field  in  the 
depletion  region  is  along  the  y  axis  and  in  the  channel  region  is  along  the  z  axis  [7]. 

The  depletion  width  h(z)  at  any  region  along  the  channel  is  infered  from  equation 
2  38' 

where  V (z)  is  the  channel  potential  at  point  z  arising  from  the  finite  drain  bias.  So 
at  the  source  end,  we  have: 


h(z  =  0)  = 


2  e(VG  +  Vu)' 
qNo 


1 

2 


2e(Va  +  Va  +  V^j 
qND 


1 

2 


(2.40) 


(2.41) 


h(z  =  L)  = 
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At  h(z=L)=a,  the  channel  is  pinched  off  and  saturation  of  the  current  occurs.  There¬ 
fore,  the  pinch-off  voltage  is  given  by: 


qNna1 


(2.42) 


This  is  the  voltage  needed  between  gate  and  drain  to  result  in  saturation  of  the 
current.  It  is  also  the  voltage  between  gate  and  source  to  shut  off  the  current  since 
we  achieve  complete  channel  depletion  at  that  voltage.  The  threshold  voltage  is 
therefore  defined  as: 

Vth  =  VP-  Vu  (2.43) 

Since  the  field  is  small  in  this  approximation,  there  is  only  a  small  deviation  to 
quasi-  neutrality  and  thus  only  small  gradient  in  carrier  concentration.  To  the  first- 
order  approximation,  we  may  neglect  any  diffusion  current.  So: 


Id  =  qNDnW[a  -  h(z)]- 


(2.44) 


where  W  is  the  width  of  the  device.  By  squaring  the  equation  for  h(z)  and  differen¬ 
tiating  it,  we  obtain: 

2hdh  =  ~dV  (2.45) 

qND 


and  therefore 


Id  =  WqNDp{a  -  h)^-h~ 

£  az 


(2.46) 


By  applying  the  current  continuity  condition  and  integrating  over  the  above  equation, 
we  obtain; 

,2.47) 

Id  =  -  *2(0)1  -  |(A3(i)  -  A3(0)j}  (2.48) 

Figure  2.5  shows  an  ideal  I-V  characteristic  for  the  MESFET.  The  current  is  saturated 
when  the  gate  to  source  voltage  reaches  Vp. 
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Figure  2.5:  Ideal  I-V  Characteristics  of  MESFET. 


2.5.3  Numerical  Solution  of  MESFET 

In  order  to  obtain  an  accurate  solution  for  the  electrical  characteristics  of  the  MES¬ 
FET  that  is  applicable  to  short  channel  devices,  which  are  most  often  used  in  actual 
applications,  we  must  solve  the  set  of  basic  semiconductor  device  equations  numeri¬ 
cally  along  with  suitable  boundary  conditions. 

Let’s  write  the  potential  $  in  the  Poisson’s  equation  as  the  intrinsic  Fermi  level. 
Then  [5]: 


(2.49) 


As  a  result,  the  boundary  conditions  at  the  source  and  drain  region,  which  are 
ohmic  are: 

"  =  5(W  -  K)  +  \/W  -  +  4"?  (2.50) 


T  ..  kT  n 
*  =  Va  +  —  In  — 
q  Hi 


where  Va  is  the  applied  voltage. 


(2.51) 
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For  the  Schottky  gate  contacts,  we  have  the  following  relations: 


Eg  kT  Nc  .. 

®  =  X  +  Tp  +  "S-  In— —  <f>m  +  Va 
2  q  2  q  N„ 


(2.52) 


(  JC*  +  Ant  +  C\ 
n  =  -qvn \n - - - J 


(2.53) 


where  <f>m  is  the  metal  work  fimction,  x  the  electron  affinity  of  the  semiconductor  and 
vn  the  surface  recombination  velocity. 

For  regions  interfacing  with  insulators,  the  potential  and  electron  concentration 
must  obey  the  following  relationships  at  the  boundary. 


n\u  — 

—  =0  •  n  =  0 

on 


(2.54) 


These  boundary  conditions  plus  the  basic  set  of  semiconductor  device  equations 
are  implemented  in  Pisces  [8],  a  device  simulation  program.  This  program  is  further 
refined  in  a  commercial  program  called  Medici  [9].  This  is  the  program  that  we  shall 
use  to  model  our  MESFETs. 


Chapter  3 

The  Properties  of  EL2 


The  native  defect  EL2  in  semi-insulating  GaAs  crystals  grown  by  the  liquid  encapsu¬ 
lated  Czochralski  (LEC)  method  has  been  the  subject  of  intensive  research  in  recent 
years.  This  defect,  which  is  associated  with  the  As  antisite  (Asgo)  and  is  responsible 
for  semi-  insulating  behavior,  also  plays  an  important  role  in  affecting  GaAs  MES- 
FET  device  parameters.  As  we  will  see  later,  any  variation  in  the  concentration  of 
EL2  will  effect  the  uniformity  of  electrical  parameters  of  devices  fabricated  by  direct 
ion  implantation  into  semi-inslating  substrates.  This  in  turn  will  effect  the  realiza¬ 
tion  of  high  performance  GaAs  LSI  circuits.  Despite  an  intensive  research  effort  to 
understand  the  nature  of  EL2,  much  controversy  still  remains. 


3.1  Physical  Property 

3.1.1  The  Donor  Nature  of  EL2 

The  donor  nature  of  EL2  was  clearly  demonstrated  by  Mircea  et  al.  [10]  Their  ex¬ 
periment  consisted  of  studying  EL2  outdiffusion  after  a  600-750° C  annealing  in  VPE 
materials.  The  experiment  clearly  established  that  EL2  centers  are  electrically  neu¬ 
tral  when  occupied  by  electrons  and  are  positively  charged  when  these  electrons  are 
released. 

From  the  thermal  electron  emission  rate  and  electron  capture  cross  section  one 
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can  obtain  the  free  energy  of  ionization,  Er  with  respect  to  the  conduction  band 
minimum  as:  [11] 

Et  =  (0.730  ±  0.02)  -  -  ■  (3.1) 

v  ’  T  +  204  v  ’ 

where  a  —  (3.5  ±  0.6)  x  10-4eVK-1. 

3.1.2  EL2  as  Asca 

EL2  is  thought  to  be  a  native  defect  which  contains  an  arsenic  antisite  (Asc0)  at 
least  as  part  of  its  physical  structure.  This  has  been  established  by  stoichiometry 
studies  [12]  and  electron  paramagnetic  resonance  data  [13].  To  establish  the  connec¬ 
tion  between  EL2  and  AsGa,  one  can  compare  the  spectroscopic  data  for  EL2  and 
Asgu •  The  comparison  is  complicated  by  the  fact  that  EL2  was  originally  detected 
and  defined  by  electrical  techniques,  while  -4soa  was  studied  by  magnetic  resonance. 
Fortunately,  when  suitably  combined  with  optical  methods,  both  techniques  provide 
information  about  the  optical  properties  of  the  respective  species.  This  optical  con¬ 
nection  is  considered  to  be  the  key  to  associating  EL2  with  Asg0. 

3.1.3  Metestable  State  of  EL2 

The  EL2  defect  posesses  both  normal  and  metastable  states  (EL2*).  The  transforma¬ 
tion  from  normal  to  metastable  states  can  be  achieved  by  illumination  of  l.leV  energy 
monochromatic  light  at  low  temperature.  The  metastable  state  can  be  brought  back 
to  the  normal  state  by  heating  the  sample  above  120K.  The  EL2*  — ►  EL2  trans¬ 
formation  is  believed  to  occur  through  the  neutral  state  EL2°  with  a  much  higher 
probability  than  the  charged  state  (EL2)+  [14].  The  configuration  of  EL2*  remains 
unknown  since  no  electrical  or  optical  properties  of  this  state  has  yet  been  discovered. 
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3.2  Theoretical  Models  of  EL2 

3.2.1  The  isolated  Asca  model 

The  isolated  Asaa  model  is  the  simplest  physical  model  since  the  correlation  of  EL2 
with  >4s<3a-  The  main  concern  of  supporters  of  this  theory  has  been  the  search  for 
large  lattice  relaxations.  Since  Bachelet  and  Scheffler  [15]  proposed  that  relaxations 
associated  with  the  isolated  Asa*  defect  cannot  introduce  a  LLR,  which  is  believed 
necessary  to  account  for  the  metastable  state  of  EL2. 

It  was  finally  shown  that  metastability  could  also  arise  from  the  isolated  AsGa 
defect  when  it  underwent  a  symmetry  breaking  distortion.  By  using  the  self-consistent 
Green’s  function  method,  the  metastable  state  is  called  the  M  configuration  and 
identified  as  a  gallium  vacancy  plus  an  arsenic  interstitial  (Vbo  —  As{)  pair  defect. 
However,  many  recent  experiments,  such  as  neutron  irradiation  show  that  EL2  is  a 
more  complex  defect  than  a  simple  point  defect  [14]. 

3.2.2  Arsenic  Antisite- Arsenic  Vacancy  —  Vas) 

This  model  was  theoretically  studied  by  a  Green’s  function  method.  It  was  found 
that  Voa  is  a  multiple  acceptor  which  can  be  either  neutral  or  negatively  charged. 
The  Asgo  ~  Vas  is  formed  by  the  probable  atomic  motion  occuring  during  annealing 
of  a  crystal,  which  causes  one  of  the  four  nearest  neighbor  arsenic  atoms  to  fill  the 
gallium  vacancy. 

This  model  also  satisfies  the  metastablity  property  of  EL2,  but  conflicts  with 
the  electron  paramagnetic  resonance  data  (EPR).  However,  EPR  measurements  were 
performed  only  on  electron  irradiated  materials  and  a  definitive  conclusion  requires 
additional  experiments. 

3.2.3  Arsenic  Antisite-Divacancy  ( Asca  —  Vba  —  Vas ) 

According  to  this  model,  the  EL2  configuration  is  formed  by  a  Vfe0  —  Va,  on  second  and 
third  neighbor  positions  with  respect  to  Asg„.  The  EL2*  configuration  corresponds 
to  an  Asgo  separating  the  two  vacancies. 
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Detailed  calculations  have  not  as  yet  been  made  for  this  model;  however,  the  latest 
interpretation  of  the  EL2  EPR  spectrum  supports  the  Asgq  -  Voa  -  Va»  model.  In 
addition,  recent  thermal  annealing  measurements  also  support  this  three  point  defect 
model. 


3.2.4  Arsenic  Antisite- Arsenic  Interstitial  ( A$ca  —  Asi) 

In  this  model,  the  stable  and  metastable  atomic  configurations  of  EL2  are  explained 
as  As,-  in  the  second  and  first  nearest  neighbor  positions,  respectively.  Self-consistent 
Green’s  function  techniques  are  carried  out  to  calculate  the  electronic  structure  and 
binding  energy  of  the  model.  It  is  found  that  in  order  to  identify  EL2  with  this 
model,  the  Ast  atom  has  to  be  positively  charged,  that  is  the  defect  state  has  to  be 
j4sGa  —  Asf.  This  finding  requires  that  an  acceptor  level  be  present  to  account  for  the 
compensation  mechanism.  This  acceptor  level  has  not  as  yet  been  identified.  However, 
the  AsGa  -  As,  model  is  supported  by  EPR  and  deep-level  transient  spectroscopy 
data.  In  addition,  optically  detected  electron  nuclear  double  resonance  (ODENDOR) 
experiments  on  semi-insulating  GaAs  material  also  finds  a  second  defect  near  Asoa, 
which  is  identified  as  the  arsenic  interstitial  [14]. 

3.3  The  Metastable  State  EL2* 

The  understanding  of  EL2*  is  one  of  the  keys  to  understanding  the  nature  of  EL2, 
unfortunately,  when  the  stable  state  is  transformed  into  the  metastable  state,  its 
normal  optical  and  electrical  properties  axe  quenched  out  and  the  state  becomes 
virtually  inaccessible  experimentally.  As  a  result,  very  little  is  known  about  E12*. 

Since  no  electrical  or  optical  properties  are  known  about  E12*,  studies  of  E12’  are 
concentrated  on  the  recovery  to  the  stable  state,  E12m  — »  EL2.  It  is  known  that 
the  recovery  of  E12*  can  be  achieved  by  heating  the  sample  to  about  120K  or  higher 
for  a  short  period  of  time.  The  activation  energy  for  this  thermal  process  has  been 
measured  by  various  workers  and  their  values  vary  from  about  0.25  to  0.36  eV.  This 
indicates  the  possibility  of  sample  dependent  activation  energy,  which  is  confirmed  by 


CHAPTER  3.  THE  PROPERTIES  OF  EL2 


27 


further  experiments.  This  indicates  that  the  metastable  state  does  not  have  a  unique 
structure  and  lends  support  to  the  concept  of  the  “EL2  family”  of  defects. 

The  lack  of  any  known  electrical  or  optical  properties  has  led  to  the  proposal 
that  E12*  is  electrically  and  optically  inactive.  A  direct  consequence  of  non-electrical 
activity  is  that  EL2  dominated  semi-insulating  GaAs  should  become  p- typed  after 
EL2  is  transformed  into  its  metastable  state. 

As  we  know  that  the  semi-insulating  behavior  of  GaAs  is  produced  by  compen¬ 
sation  of  a  residual  shallow  acceptor  by  the  deep  level  EL2,  effectively  pinning  the 
Fermi  level  near  mid  gap  which  produces  high  resistivity.  If  EL2  is  rendered  electri¬ 
cally  inactive,  then  the  electrical  properties  of  the  material  should  be  dominated  by 
the  shallow  acceptors.  This  means  that  a  thermal  current  should  be  detectable  for 
all  temperatures  above  8K  until  EL2  recovers  from  the  metastable  state. 

Experiments  were  performed  by  Jimenez  [16]  in  accordance  with  this  theory.  He 
observed  a  persistent  photocurrent  (PPC).  However  the  author  stated  that  the  PPC 
effect  was  not  universally  observed  and  only  observed  for  certain  samples.  In  addition, 
the  reported  hole  concentration  of  4  x  109  cm-3  is  too  low  to  be  thermally  ionized 
carbon  at  their  experimental  temperature  of  77K. 

A  survey  of  remainder  of  the  literature  on  this  subject  indicates  that  the  occurence 
of  the  PPC  is  an  exception  rather  than  the  rule  as  required  by  the  non-electrical 
activity  hypothesis  of  EL2.  The  question  naturally  arises  that  if  there  is  not  any 
PPC  after  quenching  of  EL2,  the  acceptors  must  be  compensated  by  a  donor  state. 
The  simplest  answer  is  to  assume  that  E12*  is  responsible  for  such  compensation. 
Calculations  of  electronic  properties  of  the  Asoa  —  As,  model  and  isolated  Asaa 
model  shows  such  possiblities  do  exist;  however,  further  calculations  are  needed  to 
give  a  difinitive  answer. 
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3.4  Optical  Mapping  of  EL2 

3.4.1  Near-infrared  Transmittance  Mapping 

For  a  wafer  of  uniform  thickness  t  with  polished  planar  faces,  the  reflection  coefficient 
for  a  paraxial  incident  flux  of  light,  using  focusing  lens  of  large  f  number  is: 

« 

where  n  is  the  refractive  index  and  is  fitted  to  experimental  data  at  room  temperature 
as:  [17] 

naooihv)  =  {7  10  +  3.78[1  -  0.180(hi/)2]-1}1/2  (3.3) 

Let  the  paraxial  incident  photon  flux  be  Io  and  flux  emerging  from  the  far  side  of  the 
wafer  be  It,  then  the  fractional  transmittance  can  be  shown  to  be: 

T  =  !f  =  (l-R)2  exp(-crf)[l  -  R 2 exp(— 2 of)]*1  (3.4) 

I o 

where  a  is  the  optical  absorption  coefficient  and  can  be  solved  for  as: 

exp(ot)  =  +  II  +  4R2T2(l  -  flnl/2}  (3.5) 


Using  the  above  equations,  for  a  sample  of  known  thickness  t  and  a  selected  hi/,  the 
measured  value  of  T  will  immediately  yield  a  value  for  a. 

For  undoped  GaAs  containing  EL2  at  a  concentration  N°  in  the  neutral  state  and 
7V+  in  the  ionized  state,  the  total  absorption  coefficient  for  near  IR  extrinsic  hu  can 
be  expressed  as: 

*  =  a  °nN°  +  <r°pN+  +  o0  (3.6) 

where  <7°  is  the  optical  absorption  cross  section  for  neutral  EL2  sites,  measured  by 
photocapacitance  methods  by  Chantre  et  al  [18],  is  the  hole  photoionization  cross 
sections,  measured  by  P.  Silverberg,  et  al  [19].  The  final  term  qq  lumps  together  any 
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absorption  processes  unrelated  to  EL2,  which  must  be  quite  small  at  room  temper¬ 
ature,  at  which  the  experiments  are  performed.  EL2  mapping  experiments  such  as 
this  has  been  conducted  by  P.  Dobrilla,  et  al.  [17,  20]. 

3.4.2  Mapping  by  Photoluminescence 

The  results  obtained  using  the  near-IR  transmittance  method  have  a  limited  spa¬ 
tial  resolution  on  the  order  of  the  thickness  of  the  wafer.  Furthermore,  they  are 
the  average  value  over  the  thickness  of  the  wafer,  and  may  or  may  not  reflect  the 
surface  device  active  area  of  the  wafer.  This  shortcoming  is  overcome  by  using  the 
Photoluminescence  (PL)  method. 

When  light  of  photon  energy  larger  than  the  band  gap  is  incident  on  the  surface  of 
the  wafer,  excess  carriers  are  created  in  the  conduction  and  valence  bands.  Carriers 
of  one  type  can  recombine  directly  with  those  of  the  other  type,  or  be  trapped  first  by 
a  deep  level  and  then  recombine.  PL  occurs  as  a  result  of  the  radiative  recombination 
process. 

Let’s  assume  that  the  carrier  lifetime  (t)  is  limited  by  recombination  via  “lifetime- 
killer”  centers,  and  the  concentration  of  excited  carriers  (nex)  is  larger  than  the  intrin¬ 
sic  carrier  concentration  and  the  concentration  of  the  deep  levels  Ng.  The  intensity 
of  the  band-to-band  transition  (/g)  and  that  of  the  deep  level  (Id)  are  [21]: 

/g  a  n2exr  (3.7) 

Id  °c  nexrND  (3.8) 

Therefore,  under  a  constant  excitation  condition,  /g  is  proportional  to  r,  or  inversely 
proportional  to  the  concentration  of  the  “killer  centers.”  This  implied  that  the  spatial 
variation  of  7g  reflects  the  spatial  inhomogeneity  of  the  crystalline  quality.  The 
intensity  Id ,  on  the  other  hand,  is  proportional  to  tNd ■  Assuming  that  the  variation 
of  r  on  a  wafer  is  small,  then  the  variation  of  Id  reflects  the  distribution  of  the  deep 
levels. 

Such  PL  experiments  were  conducted  by  Tajima  [22].  The  measurement  was 
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carried  out  at  room  temperature  in  air.  The  sample  was  mounted  on  an  X-Y  stage 
with  a  travel  distance  of  100  mm  for  both  axes.  The  diameter  of  the  laser  beam 
varied  from  0.1  to  1  mm  with  an  incident  power  of  0.1  to  100  W  cm-2.  The  PL 
signal  is  analyzed  with  an  apparatus  of  high  optical  efficiency,  consisting  of  a  grating 
monochromator  with  a  75  and  a  600  grooves/mm  grating  blazed  at  2.0  and  1.0  fim, 
respectively,  and  Ge  and  PbS  detectors.  A  deep  level  PL  band  with  a  peak  of  0.65 
eV  was  observed  and  is  believed  to  be  traceable  to  EL2. 


3.5  Effect  of  EL2  on  Device  Parameters 

The  effect  of  dislocations  and  EL2  on  the  electrical  properties  of  devices  have  been 
widely  studied  and  debated.  Although  EL2  as  a  donor  is  electrically  inactive  in  the 
n  channel  FET  fabricated  by  direct  ion  implantation,  it  can  well  affect  the  activation 
ratio  of  the  implants,  especially  the  amphoteric  Si  implant,  thus  affecting  the  device 
parameters.  Studies  performed  by  Dobrilla  and  Blakemore  [23]  show  that  a  large  value 
of  N°,  the  neutral  EL2  concentration,  tends  to  make  Id.,  larger  and  Vth  more  negative. 
Their  experiments  are  done  using  near-IR  transmittance  mapping  of  EL2  and  KOH 
etching  to  reveal  dislocation  etch  pits  following  various  electrical  measurements.  They 
report  that  one  of  their  wafers  had  a  uniform  dislocation  density  as  revealed  by  KOH 
etching  and  a  radially  symmetrical  N°  distribution.  They  were  able  to  correlate  the 
device  parameters  with  EL2  density  for  this  particular  wafer. 

However,  studies  done  by  Miyazawa  and  Wada  [24]  found  evidence  for  the  influence 
of  dislocation  proximity  to  a  device  on  threshold  voltage.  They  found  an  increase  of 
5  ±  2 mV  in  Vth  per  1/rm,  for  a  gate  to  pit  distance  below  about  50jim.  It  was  also 
found  that  E12  and  dislocations  are  correlated  [25].  The  results  by  Miyazawa  and 
Blakemore  do  not  in  fact  contradict  each  other. 


Chapter  4 

Experimental  Setup 


We  have  taken  an  extremely  versatile  approach  to  focus  on  the  study  of  issues  effecting 
GaAs  MESFET  electrical  characteristics.  Our  purpose  is  to  develop  and  substantiate, 
through  experimental  design,  fabrication  and  testing,  models  for  understanding  of  the 
relationships  between  material  and  process  defects,  or  variations,  on  the  performance 
and  uniformity  of  GaAs  digital  ICs.  Our  approach  is  to  couple  the  modeling  and 
design  effort  with  an  advanced  manufacturing  facility  to  allow  controlled,  state-of- 
art  device  fabrication  processes  to  be  used  to  test  the  models  being  deveploped  and 
to  assist  in  the  design  of  realistic  test  structures.  Therefore,  a  unique  part  of  this 
effort  is  the  cooperative  design  of  test  structures  specifically  aimed  at  substrate  ma¬ 
terial  defects  and  ion  implantation  pr  cess  variations  using  a  state-of-art  sub-micron 
MESFET  foundry  process.  Further  test  structures,  aimed  at  such  important  topics 
as  device  scaling  and  gate  orientation,  are  also  included.  These  test  structures,  are 
part  of  a  reticule  which  is  stepped  over  the  entire  wafer  to  minimize  the  effects  of 
lithographic  misalignment.  Computer  automated  testing  of  the  test  structures  allows 
extension  of  the  modeling  effort  to  wafer  scale,  which  due  to  the  inherent  foundry 
process  control,  allows  the  eventual  prediction  of  wafer  yields  resulting  from  spatially 
dependent  defects  and  processes,  such  as  dislocations  and  ion  implantation. 
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Table  4.1:  Major  Structures  and  Design  Rules 


Major  Structures 

Design  Rule 

FET  Array 

Size  l/i  x  l/i  ,  gate  to  gate  spacing  11.5  fxm. 
E-D  cell  spacing  22.5  /im 

D-mode,  E-mode 

Gate  length:  0.5,  1.2,  2.0,  5.0,  8.0  ftm. 
Width:  2,  10,  20  and  50  /im 

4.1  Dense  Array  Pattern  Structures 

The  Advanced  Test  Structure  mask  set  was  designed  for  use  with  the  state-of-the-art, 
self-aligned,  refractory  gate  MESFET  process  available  at  the  Rockwell  foundry.  The 
test  structures  utilize  fabrication  technology  capable  of  producing  both  enhancement 
(E)  and  depletion  (D)  mode  MESFET  devices  on  the  same  wafer.  We  however  utilize 
this  technology  to  implant  different  ion  species  into  the  channel  of  MESFETs  in  close 
proximity.  This  allows  a  direct  comparison  of  the  effects  of  different  ion  species  under 
similiar  conditions.  This  approach  eliminates  complications  caused  by  spatially  or 
wafer  dependent  defect  structures  or  background  acceptor  concentrations. 

Table  4.1  identifies  the  various  structures  available  on  the  ATS  for  testing.  As  seen 
in  this  table,  the  ATS  contains  channel  implants  of  Si  and  Se  dopants  in  FETs  in 
close  proximity.  FETs  with  different  gate  lengths  and  gate  widths  are  also  fabricated 
for  examination  of  short  channel  and  backgating  effects.  Gate  lengths  include  0.5, 
0.6,  0.7,  0.8,  0.9,  1.0,  1.1,  1.2,  2,  5  and  8  /im. 

4.2  Wafer  Processing  Conditions 

The  following  list  outlines  the  fabrication  process  for  what  we  shall  refer  to  as  Lot 
1  wafers.  There  are  a  total  of  fourteen  (100)  oriented  wafers  taken  from  the  same 
boule,  processed  in  this  Lot. 

•  n~  channel  implants. 


1.  Se:  Energy=124  KeV,  Dose=3.33El2  cm~2. 
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2.  Si31:  Energy=50  KeV,  Dose=4.75El2  cm-2. 

3.  SiF :  Energy =85  KeV,  Dose=4.75El2  cm-2. 

•  Anneal:  Temperature=850°C,  Time=20  minutes. 

•  Gate  Lithiography. 

•  LDD  nitride  deposition. 

•  n+  implant  of  SiF,  Energy=150  KeV,  Dose=2.2El3  cm-2. 

•  Anneal:  Temperature=790°C,  Time=20  minutes. 

•  isolation  implant. 

•  ohmic  lithiography. 

•  alloying. 

•  nitride  deposition. 

•  via  definition. 

•  metal  1. 

Figures  4.1  and  4.2  show  the  dopant  profiles  of  Lot  1  process.  These  profiles  are 
obtained  from  the  Stanford  GaAs  process  simulator,  Suprem3.5  [37]. 

For  Lot  2,  we  repeat  the  Lot  1  process  except  that  wafers  from  three  different 
vendors  are  used.  Since  the  absolute  value  of  the  threshold  voltage  for  the  Se  channel 
implanted  devices  in  Lot  1  was  found  experimentally  to  be  below  target,  we  increased 
the  Se  implant  dose  from  3.33  x  1012cm-2  to  4.00  x  1012cm-2  for  this  Lot. 

Figure  4.3  and  4.4  show  the  threshold  voltage  maps  across  two  wafers.  We  notice 
that  there  is  quite  a  bit  of  threshold  voltage  variation  across  the  wafers. 
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Lot  1  Si  Implant  Profile 


Figure  4.1:  Si  profiles  calculated  for  the  Lot  1  process  using  Suprem3.5. 


Lot  1  Se  Implant  Profile 


Figure  4.2:  Se  profiles  calculated  for  the  Lot  1  process  from  Suprem3.5. 
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Figure  4.4:  Threshold  voltage  mapping  for  a  Lot  1  Se  implanted  wafer. 
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Figure  4.5:  Layout  for  each  3  inch  wafer. 


4.3  Wafer  Testing 

Each  3  inch  wafer  is  divided  into  44  usable  grids  and  identical  devices  are  fabricated 
in  each  grid.  The  wafer  layout  in  shown  in  Figure  4.5. 

The  devices  selected  for  testing  including  gate  lengthes  of  0.5,  0.7,  1.1,  2.0,  5.0/un 
and  device  oriented  at  0,  ±45,  90  degree  to  the  major  axis  for  both  implant  species. 
This  totals  to  40  distinct  devices  tested  in  each  grid.  Half  of  the  44  grids  are  stepped 
over  via  an  automatic  testing  probe.  As  a  result,  a  total  of  880  MESFETs  are  tested 
on  each  wafer  and  the  resulting  threshold  voltages  are  compiled  in  formatted  ascii 
files.  For  each  lot,  threshold  voltage  data  is  obtained  for  12320  descrete  devices. 


Chapter  5 

Equilibrium  Properties  of  Pure 
GaAs 


5.1  The  Thermodynamics  of  Chemical  Reactions 


Under  constant  pressure,  equilibrium  is  characterized  by  a  minimum  in  the  Gibbs 
free  energy.  Let  us  consider  the  following  chemical  reaction: 


<M(s)  +  bB{s)  ~  cC(s)  +  dD(s )  AG? 
The  change  in  Gibbs  free  energy  is  of  course  given  by: 

AG  =  AH  -  TAS 


(5.1) 


(5.2) 
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Let  p°A,  p°B,  p^  and  p°D  be  the  equilibrium  partial  pressures  of  the  pure  components. 
Then  we  may  convert  the  above  reaction  (5.1)  into  one  involving  the  gaseous  compo¬ 
nents  of  the  constituents  only,  [26]. 

aA(s)  +  f>B(s)  =  cC(a)  +  dD(a)  A  G? 

aA(g,  pA  )  =  aA(s)  AGj  =  0 

aB(g,p°B)  =  aB(s )  AG3  =  0  (5.3) 

cC(s)  =  cC(g,p°c)  AG  4  —  0 

dD(s)  =  dD(g,p°D)  AG2  =  0 

Adding  up  the  above  equations,  we  obtain: 

aA(g,p° )  +  bB(g,p°B)  =  cC(g,p°c )  +  dD(g,p°D)  A G6  =  AG?  (5.4) 

The  free  energy  change  A G2  to  AGs  are  all  zero  becasue  the  p°  are  chosen  to  be 
the  equilibrium  values  for  the  pure  components.  Let  us  now  choose  another  set  of 
partial  pressures  peA,Ps,  Pc >  Pd  so  that  the  reaction  5.4  will  be  in  equilibrium.  Hence 
the  change  in  Gibbs  free  energy  will  be  zero. 

aA(9,PA)  +  bB(g,p%)  =  cC(g,pec )  4-  dD(g,peD)  A G7  =  0  (5.5) 

Next,  we  need  to  calculate  the  change  of  Gibbs  free  energy  for  each  pure  substance 
from  a  partial  pressure  of  p°  to  p‘  so  that  we  can  find  a  relationship  between  AG?  to 
A Gr.  The  thermodynamic  function  satisfied  by  the  Gibbs  free  energy  is: 

dG  =  Vdp  -  SdT  (5.6) 


Let  the  change  of  partial  pressure  take  place  under  constant  temperature  and  further 
assume  the  vapors  to  be  an  ideal  gas,  then  pV  =  RT  and  we  obtain: 

RT 

dG  =  Vdp  =  —dp 
P 


(5.7) 
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AG  =  RTln^r 
P? 

Substituting  A,  B,  C  and  D  into  the  above  equation,  we  obtain; 


aA(9iP°A)  +  bB(g,p%)  =  cC(g,p%)  +  dD(g,p°D) 
«A(9,Pa)  =  aA(g,p°A) 

*>B(g}p%)  =  bB{g,p%) 

<£(9,  Pc)  ~  cC(9,Pc) 
dD(g,p°D)  =  dD(g,peD) 

aA^PeA)  +  bB(g,p'B)  =  c  C(g,p'c)  +  dD(g,p‘D) 
Upon  summation  of  all  the  free  energy  changes, 


AGe  =  AG° 
AG8  =  aRT\n& 

AG»  =  bRT  In  4 

p% 

AGio  =  c/Crinj| 

AGn  =dRT\n% 

PD 

AGr  =  0 


AG7  =  0  =  AG6  +  AG8  +  AGg  +  AGi0  +  AGU 


0  =  A  G°  +  RT  In  ^cIPcYiPblPpY 
(^M)°(Pb/Pb)6 


(5.10) 


(5.11) 


We  see  that  the  argument  of  the  logarithmic  function  depends  only  on  the  temper¬ 
ature  of  the  reaction  and  is  independent  of  a,  b,  c,  and  d.  We  call  this  the  equilibrium 
constant  K. 


AG°  =  -RT\nK 

Let’s  define  the  activity  of  a  component  as: 


(5.12) 


„  -  Pi 


Then,  the  equilibrium  constant  can  be  written  as: 

K  =  (ac,)c(ap)‘i  _  n(q,?)*n,-)(Products) 
(aA)a(aB)b  n(a*)^n,)(Reactants) 


(5.13) 


(5.14) 
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5.2  Dilute  Solution  and  Henry’s  Law 

Consider  a  solution  of  n  solute  atoms  and  N-n  solvent  atoms.  Define  the  atom  fraction 
as: 

X  =  Jr  (5.15) 

Let  the  partial  vapor  pressure  of  solute  atoms  be  p,-  and  the  vapor  pressure  of  the 
pure  substance  be  pf.  Then  if  X  =  1,  obviously  we  would  have  p,  =  p°.  On  the  other 
hand,  if  X  =  0,  then  p,  =  0.  An  ideal  solution  is  defined  as  one  in  which  the  vapor 
pressure  of  a  component  i  is  a  linear  function  of  composition: 

Pi  =  XiP ?  (5.16) 

This  is  known  as  Raoul t’s  Law.  Comparing  this  with  equation  5.13: 

=*  <ii  =  Xi  (5.17) 

Most  real  solutions  do  not  fall  into  this  catatory.  However,  as  long  as  the  separation 
between  solute  atoms  are  large  and  the  interaction  is  small,  we  may  use  the  dilute 
solution  approximation. 

Pi  =  bXi  (5.18) 

The  above  equation  suggests  that  the  thermodynamic  properties  of  the  solute  are 
additive  for  dilute  solutions.  Thus,  the  vapor  pressure  of  a  solution  containing  two 
hundred  atoms  would  be  twice  that  containing  one  hundred  atoms.  Here  b  is  a 
constant  whose  value  depends  on  the  nature  of  the  solute  and  solvent.  This  is  known 
as  Henry’s  Law. 

Again,  from  the  definition  of  activity,  we  obtain: 

ai  =  -jgXi  or  a{  =  7?X,  (5.19) 

where  7  is  called  the  activity  coefficient,  7,  =  1  for  an  ideal  solution.  It  is  constant 
for  a  given  solute  in  a  dilute  solution  independent  of  concentration. 
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For  an  dilute  solution,  the  equilibrium  constant  in  equation  5.14  is  given  by: 

,,  _  (XcnxpY  _  \CY\Pl 

where  the  square  brackets  indicates  the  atomic  concentration  of  an  element  and  the 
activity  coefficient  has  been  absorbed  into  the  equilibrium  constants. 


5.3  Defect  Nomenclature 

Three  basic  types  of  defects  may  occur  in  a  crystal.  An  unoccupied  lattice  site  is 
called  a  vacancy.  An  atom  occuping  a  non-lattice  site  is  called  an  interstitial.  In 
addition,  foreign  atoms  can  occupy  either  the  vacancy  or  interstitial  sites.  Table  5.3 
lists  the  defect  nomenclature. 

When  writing  reaction  equations  involving  defects,  three  considerations  must  ap¬ 
ply  as  follows  [26]: 

1.  Total  charge  must  be  conserved  on  both  sides  of  the  equation. 

2.  Mass  must  be  conserved  on  both  sides  of  the  equation. 


CHAPTER  5.  EQUILIBRIUM  PROPERTIES  OF  PURE  GAAS 


43 


3.  Fixed  proportion  of  Ga  and  As  sites  must  be  created  regardless  of  whether  or 
not  they  are  actually  occupied  or  not. 

For  a  compound  semiconductor  such  as  GaAs,  the  site  ratio  between  Ga  and  As  as 
well  as  the  gross  composition  must  be  preserved;  therefore,  different  types  of  defects 
mn  not  occur  arbitrarily.  There  are  six  basic  types  of  combinations  as  follows  [27]: 

1.  Schottky  disorder  involving  vacancies,  Vfea  and  V 

2.  interstitial  disorder  involving  interstitial  atoms,  Ga;  and  As;. 

3.  anti-structure  disorder  involving  misplaced  atoms,  Ga a,  and  Asc0- 

4.  Schottky  interstitial  disorder  involving  vacancies  and  interstitials  for  the  same 
type  of  atom,  e.g.,  Voa  +  Ga,  or  Va*  +  As,.  This  type  of  defect  combination  is 
usually  called  a  Frenkel  defect. 

5.  Schottky  anti-Site  disorder  involving  vacancies  and  misplaced  atoms  for  the 
same  type  of  atom,  e.g.,  Vg„  +  Goa,  or  Va,  +  Asg«. 

6.  Interstitial  anti-site  disorder  involving  interstitial  atoms  of  one  type  and  mis¬ 
placed  atoms  of  the  other  type,  e.g.,  Ga;  -f  Asc0  or  As;  +  Ga^,. 

5.4  Native  Defects  of  GaAs 

5.4.1  Schottky  Defects 

For  semi-insulating  GaAs,  there  are  usually  more  As  than  Ga  atoms.  If  Schottky 
defects  are  the  only  disorder,  then  there  must  be  an  equivalent  number  of  Ga  vacancies 
to  accomodate  for  this. 

Nv  =  Na,  —  Ng<i  (5.21) 

In  thermal  equilibrium,  vacancies  of  both  type  will  be  formed  due  to  the  entropy 
change  brought  about  by  the  vacancies. 


VGa  +  Va,  ^  0 


(5.22) 
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Assume  Mv  vacancies  of  Ga  and  As  are  formed  through  this  process,  the  generaJ 
form  of  the  Gibbs  free  energy  is  given  by  [27]: 

G  =  NoaM(Ga)  +  NAtM(As)  +  (Ny  +  My)M(Ga)  +  My/i(As)  -  TSconJt3  (5.23) 


Where  p  is  the  chemical  potential  of  the  species.  S^fig  =  k  In  W  where  W  is  the 
statistical  distribution  of  vacancies  on  the  available  sites  and  is  the  product  of  two 
the  factors: 


1.  NGa  atoms  and  (Ny  +  Mv)  vacancies  distributed  over  (Ng0  +  Ny  +  My)  sites. 

2.  Na,  atoms  and  My  vacancies  distributed  over  (NAt  4-  My)  sites. 


Then: 

W  —  W  W  —  ^ A"  *  My )l(NGa  +  Ny  +  My )! 
1  2  “  NGa\(Ny  +  My)\NA.\My\ 

Using  Stirling’s  approximation.  In  N\  =  N  In  JV  -  N,  we  find: 


(5.24) 


Sctmjig  —  ~^(NGaln[GaGo]  +  (iVv4'Mv)ln[VGo])  —  (NAt ln[AsA*]+Afv ln[V)i»])  (5.25) 
where  the  square  brackets  indicate  the  site  fractions. 

[VGa]  = 


Ny  •+  My 

(5.26) 

N Ga  +  Ny  +  My 

Mv 

NAs  -f  My 

(5.27) 

The  number  of  vacancies  formed  in  thermal  equilibrium  is  determined  by  minimizing 
the  Gibbs  free  energy: 

(  o77-)  =M(Vba)  +  ^ln[UG0]  +  p(VrA,)  +  ^ln[V,/i,]==0  (5.28) 


Resulting  in 


[Vca  ](VAt]  =  Ks 


(5.29) 
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where 

Ks  =  exp(-(^(Voa)  +  f*(VAi))/kT)  (5.30) 

The  exact  same  relationship  can  also  be  obtained  by  mass  action  considerations. 


5.4.2  Interstitial  Defect 

The  interstitial  defect  can  be  described  by  the  following  reaction: 

Gaea  +  As a»  Gai  -f  As{  (5.31) 


According  to  the  mass  action  law: 

=  Ki  (5.32) 

Since  the  amount  of  defect  atoms  are  much  smaller  than  the  amount  of  Ga  or  As 
atoms.  The  site  fractions  of  Ga  and  As  are: 


[GaMAsi] 

[GaGa][AsAa] 


[GaGo]  ~  1  and  [As*,]  «  1 


(5.33) 


Resulting  in  the  approximation: 


[Ga,][As,]  =  Ki 


(5.34) 


5.4.3  Anti-site  Defect 

The  anti-site  defect  plays  an  important  role  in  the  electrical  properties  of  GaAs. 
The  deep  donor  EL2  is  assoicated  with  the  As  anti-site,  Asg<i,  and  is  responsible  for 
compensating  the  residual  carbon  acceptors  to  give  LEC  GaAs  its  semi- insulating 
property. 

The  equilibrium  of  anti-sites  is  described  as: 


Gaea  +  Asas  Ga, 4,  -f  Asgo 


(5.35) 
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=  Ka  (5-36) 

Again,  for  small  amount  of  defects,  [ Gaca ]  =  [Asa*]  w  1,  and  the  above  equation 
simplifies  to: 

[GaA*][AsGa\  =  Ka  (5.37) 

Once  again,  a  simple  product  relationship  exists  for  anti-site  structures  as  with 
Schottky  and  interstitial  disorder. 


[Goa*  ]  [Asca] 
[GaGa][AsA,] 


5.4.4  Frenkel  Defect 


This  type  of  defect  was  first  observed  by  Frenkel  and  arises  from  the  transfer  of  atoms 
from  lattice  to  interstitial  sites.  It’s  the  most  important  of  the  asymmetrical  type  of 

defect. 

Gogo.  ^  Gai  +  Vgo 

(5.38) 

As  a*  **  As,  4-  Va, 

(5.39) 

From  the  mass  action  law: 

(Co.Hvy 

I  Gao.]  ' 

(5.40) 

l^illKu]  ,, 
l  Asa*]  ~  2 

(5.41) 

Approximating  the  site  fraction  of  Ga  and  As  to  be  1,  results  in: 

=»  (G'a.HVca]  =  I<1  and  [ASi][VAa]  =  K2  (5.42) 

Multiplying  the  above  two  equations  together  and  using  equation  5.28,  we  obtain: 


[Co.H/M  =  ^ 


(5.43) 
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5.4.5  Schottky  Anti-site  Defect 

The  hybrid  defect  of  the  Schottky  and  anti-structure  disorder  is  described  by: 


GdOa  ^  Gaj i,  +  2  Voa 

(5.44) 

As  At  ^  ASGa  +  TV  At 

(5.45) 

[GaAt\[Vaa?  r„ 

[Gao.) 

(5.46) 

\AsgMa,V  ,, 

[ASAt]  “2 

(5.47) 

5.4.6  Interstitial  Anti-site  Defect 

The  hybrid  defect  of  the  interstitial  and  anti-site  defect  is  described  by: 


5.5  Equilibrium  Properties  of  GaAs 

5.5.1  Pure  GaAs 

Due  to  the  compund  nature  of  GaAs,  the  intrinsic  defect  structure  is  much  more 
complex  than  an  elemental  semiconductor  such  as  Si.  In  GaAs,  there  are  two  types 
of  vacancies,  two  kinds  of  interstitials  and  two  kinds  of  anti-site  structures  to  be  ac¬ 
counted  for.  The  following  reaction  describes  the  equilibrium  properties  of  crystalline 
GaAs. 


0=*e~  +  h+ 


(5.51) 


CHAPTER  5.  EQUILIBRIUM  PROPERTIES  OF  PURE  GAAS 


48 


0  Vca  +  Va.  (5.52) 

MS.  5=4  Voa  +  (5.53) 

va,  —  Y+  +  «"  (5.54) 

GaGa  +  VA,  ^  GaA.  +  Va*  (5.55) 

Asa,  +  VGa  **  AsGa  +  VA,  (5.56) 

GaGa  **  Ga  +  VGa  (5.57) 

Asa,  Asi  +  VA,  (5.58) 

\As4(g)  +  V%t  *±  As°a,  (5.59) 


The  mass  action  relationships  for  the  above  equations  are: 


ie-p+l  =  Ki 

(5.60) 

[KfcHKu]  =  Ks 

(5.61) 

(*+]lvo.l  „ 

I/O  ~ 

VGa 

(5.62) 

j, 

yo  —  n.d 

(5.63) 

[GaA,)[VG B] 

IK..]  - 

(5.64) 

(AscIIVa.  ]  „ 

[Ifcj  ~K-J 

(5.65) 

[G«()[Va]  =  AV, 

(5.66) 

lAs,]|K,.j  =  A,; 

(5.67) 

IK..]  = 

(5.68) 
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One  must  also  consider  the  charge  neutrality  condition: 


|e-]  +  |V£J  -  [A+]  +  IK£] 


(5.69) 


The  results  of  this  is  that  we  have  10  equations  and  11  unknowns.  We  will  express 
all  the  quantities  in  terms  of  the  As  pressure,  px»-  There  are  three  regions  of  interest 
to  us.  If  the  material  is  p-type  then  the  concentration  of  electrons  and  the  ionized 
donors  is  neglected.  For  semi-insulating  material  case  [e~]  =  [h+],  and  for  the  n-type 
material  case,  [V^]  =  [e— ].  These  regions  are  summarized  as  follows: 


Region  One  :  [V'J'J  =  [h+\  p  type  material 

Region  Two  :  [e~]  =  [A+]  Intrinsic  material 

Region  Three  :  [V£]  =  [e~]  n  type  material 


(5.70) 


Assume  x  =  Ka»P\u .  The  neutral  defects  can  be  solved  independent  of  the  regions 
and  are  given  as: 

IKu]  =  *_l  (5.71) 


[VGo\  =  Ksx 
\Gaju }  = 

[AsGa]  =  Ka2Ksx 2 
[G«i]  =  *-1 

[As,]  —  Ki2x 


(5.72) 

(5.73) 

(5.74) 

(5.75) 

(5.76) 


In  region  1,  the  follwoing  equations  hold: 


[A+l  =  (VS.!  =  Ji kJTsx'I* 

|£'J  -  vlk*-'1' 


(5.77) 


(5.78) 
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In  region  2,  we  have: 


Kiy/KJG 

Ki 


x-1'2 


[e'l  =  [A+]  =  JK; 

KJ  = 

™  -  vk*-' 


In  region  3,  the  solutions  are: 


(5.79) 

(5.80) 

(5.81) 

(5.82) 


(Ktl  =  [e“]  =  JZx-W 


(5.83) 


[A+]  =  —7^-X 


Ki  -1/2 

y/JT/ 


[Voa]  = 


(5.84) 


(5.85) 


If  the  natural  logrithim  of  various  concentrations  is  ploted  against  the  logrithim  of 
x,  the  graph  consists  of  straight  lines.  We  show  this  in  Figure  5.1  and  the  diagram 
is  known  as  a  Kroger- Vink  diagram  [28].  The  exact  position  of  the  lines  can  not  be 
determined  unless  we  know  the  value  of  all  the  equilibrium  constants  in  detail. 

Implanted  GaAs  is  subjected  to  high  temperature  anneals  in  which  the  equilibrium 
of  the  various  defects  are  established.  As  the  crystal  cools  down,  the  atomic  defects 
usually  don’t  have  enough  thermal  energy  to  diffuse  to  new  sites  in  order  to  establish 
new  equilibrium  at  room  temperature.  As  a  result,  the  defects  are  “frozen  in.”  The 
electronic  defects,  however  are  much  more  mobile  and  will  establish  an  equilibrium 
with  the  “frozen  in”  atomic  defects.  Therefore,  in  figure  5.1,  the  electron  and  hole 
concentrations  decrease  while  the  various  charged  defects  tend  to  become  neutral. 


5.5.2  GaAs  with  Carbon  and  EL2 

LEC  grown  semi-insulating  GaAs  contains  residual  carbon  that  is  compensated  by 
the  deep  native  donor  EL2,  that  is  associated  with  the  As  anti-site.  To  describe  the 
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logx 

Figure  5.1:  Kroger-Vink  Diagram  of  Pure  GaAs 

equilibrium  configuration  of  LEC  GaAs,  we  add  three  more  equations  to  the  general 
equilibrium  equations  discussed  for  pure  GaAs  from  the  last  section.  We  also  need  a 
new  neutrality  condition. 

\c<(g)  +  VA.  -  CAt  (5.86) 

Cm  -  C\a  +  e-  (5.87) 

Asco  —  As^a  +  h+  (5.88) 

The  neutrality  condition  is  now: 

(«']  +  [AsdJi  =  l*+l  +  [Ca.\  (5-89) 

At  elevated  temperatures  (1000K),  the  neutrality  condition  can  simply  be  approx¬ 
imated  by  [e“]  =  [/j+]  since  carbon  and  EL2  levels  are  typically  around  1  x  1016 
cm-3.  However,  at  room  temperature,  the  neutrality  condition  is  approximated  as 
[As<ja]  =  [CJJ.  In  either  case,  solutions  are  obtained  as  described  in  the  last  section. 


Chapter  6 


Activation  Model:  Si  and  Se 

Silicon  and  selenium  are  the  most  widely  used  n-type  dopants  in  GaAs.  Si  activation 
has  been  widely  studied  in  the  last  ten  or  so  years  with  Se  activation  recently  receiving 
more  attention.  Experimentally,  reported  values  for  Si  and  Se  activation  varies  widely 
from  as  low  as  40%  to  near  100%  [29,  30,  31]. 

6.1  Si  activation 

Silicon  is  an  amphoteric  dopant  in  GaAs.  It  acts  as  a  donor  when  occupying  the 
Ga  site  and  as  an  acceptor  when  occupying  the  As  site.  Therefore,  the  active  donor 
concentration  is  given  by  the  difference  between  the  Si  occupying  the  Ga  and  As  sites. 
In  our  investigation,  we  only  consider  the  charged  defects  EL2  and  carbon.  The  deep 
level  EL2  is  present  in  GaAs  in  the  range  of  around  1  x  1016  cm-3,  while  carbon  is 
typically  found  to  be  around  2  x  JO15  cm-3. 

Several  theoretical  explanations  described  the  activation  of  Si  in  GaAs.  An  early 
theory  developed  by  Greiner  is  based  on  the  following  equation  [32]. 

+  VZ  +  2e-  =  Si~As  +  Voa  +  2  h~  (6.1) 

By  utilizing  the  law  of  mass  action  and  p  =  nj/n,  the  ratio  of  Si  occupying  the  As 
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and  Ga  sites  may  be  obtained  as: 


is.;,)  nvai,«>4 

IS«SJ 


(6.2) 


Greiner  defines  the  saturation  electron  concentration,  n,at  as  the  electron  concentra¬ 
tion  occuring  when  [5*^J]/($i<sa]  =  1*  Then: 


where  Kv  is  fitted  as  a  constant  dependent  only  on  temperature. 

Although  the  charged  vacancies  can  be  neglected  in  the  neutrality  condition  since 
their  concentrations  are  typically  much  lower  than  the  intrinsic  carrier  density  at 
an  elevated  annealing  temperature.  This  does  not  imply  that  the  charged  vacancy 
concentrations  are  not  influenced  by  the  electron  concentration.  The  relationships 
between  the  ionized  and  neutral  vacancies  are: 


^  *ef-eVa.  7  (6-4) 

1  +  (j  exp(  KT**) 

i  (6-5> 

1  +  3  exp(  ■  ) 

As  the  electron  concentration  changes  due  to  different  amount  of  doping,  the  Fermi 
level  shifts  as  well.  As  a  result,  the  amount  of  ionized  vacancies  change  as  well. 
The  concentrations  of  V^a  and  change  in  opposite  directions  while  their  product 
remains  constant,  as  required  by  mass  action  law. 

In  the  theory  proposed  by  Bindal  et  al  [33],  neutral  vacancies  are  used  as  the 
starting  point  and  the  equilibrium  conditions  for  Si  on  Ga  and  As  sites  are  proposed 


Sica  +  Va»  —  Si  As  +  Voa 

[Sm.]  _  1  [Vas] 

|S*&]  K  [Van] 


(6.7) 
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Solving  this  together  with  the  total  Si  concentration,  given  as: 

[Sioa]  +  [StA.]  (6.8) 

Since  for  an  n-type  semiconductor,  most  of  the  shallow  donors  as  well  as  acceptors 
are  ionized,  the  total  Si  concentration  is  actually  given  by: 

[Si],*.,  =  [Sic]  +  [S.„.]  +  [SiJ.]  +  [s.;.]  (6.9) 

and  can  be  approximated  as: 

[St]  total  =  [Sled  +  [StA.]  (6-10) 

As  a  result,  the  mass  action  equation  proposed  by  Bindal  et  al.  is  not  complete  and 
must  also  include  terms  for  the  ionization  of  Sioa  and  Si  a,- 

Vanasupa  proposed  a  theory  based  on  the  following  mass  action  equation  [34]. 

Sit.  +  K?.  +  «-  =  Si-A.  +  VS.  +  (6.11) 


As  a  result,  the  ratio  of  Si  on  As  sites  to  those  on  Ga  sites  is  written  as: 


[Si a.)  r_  M)*9 
[StL]  '[V&]»? 


n(-r 


(6.12) 


where  fl  is  fitted  as  a  constant  dependent  only  on  temperature. 

Although  the  product  of  the  neutral  vacancies  must  remain  constant  when  equi¬ 
librium  is  reached  at  a  particular  temperature,  the  ratio  of  the  two  does  not.  This 
ratio  is  intimately  connected  with  the  stoichiometry  of  the  material  as  well  as  the 
distribution  of  dislocations. 

To  systematically  study  the  activation  of  dopants  in  GaAs,  let’s  first  list  all  the 
revelant  mass  action  equations.  The  reactions  concerning  pure  GaAs  include: 


0=±e~  +  h+ 


(6.13) 
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0  **  VCa  +  VA,  (6.14) 

\As4(g)  +  VAt  Asa$  (6.15) 

4 

The  EL2  defect  is  represented  by  Asca  and  the  revelant  reactions  include: 

As  Am  +  Vfca  5=4  AsQa  +  VA»  (6.16) 

A soa  *■*  A Sffa  +  e  (6.17) 

Equations  describing  carbon  incorporation  into  the  GaAs  lattice  are: 

C(g)  +  Vam  -  Cam  (6-18) 

Cam  -  C-At  +  h+  (6.19) 

and  the  relations  involving  Si  incorporation  into  the  GaAs  lattice  are: 

Si(g)  +  VCo  p*  SiGa  (6.20) 


Si{g)  +  V Am  5=4  Siam 
SiGa  4=4  SiQa  +  e~ 
Si  Am  *=  SiA,  +  h+ 

The  neutrality  condition  now  reads: 


(6.21) 

(6.22) 

(6.23) 


(el  +  [Cam]  +  [SCJ  =  [ h+ )  +  [AsU  +  [Si% J  (6.24) 


There  are  now  a  total  of  12  equations  and  15  unknowns.  As  a  result,  each  quantity 
may  be  expressed  as  a  function  of  3  other  variables.  Note  that  adding  each  new 
element  such  as  carbon  or  Si  adds  one  more  variable  than  equation. 

From  the  reactions  describing  the  incorporation  of  Si  into  GaAs  lattice,  we  obtain: 


[Sl'ca]  =  [VcalpSi 


(6.25) 


CHAPTER  6.  ACTIVATION  MODEL:  SI  AND  SE 


56 


[Si,.]  =  A',IKub« 

(6.26) 

[«SJ[e-]  _  ,, 

|S<£.)  L  A',[V0.)pSi 

(6.27) 

,  A'dA,1Vg.1p5. 

=>  -  [e-| 

(6.28) 

Similarily,  we  can  derive  an  equation  for  [5t'a.]> 

,  KAK,\VA.)pSi 
-  [A+] 

(6.29) 

Taking  the  ratio  of  the  above  two  equations,  results  in: 

[Si~Aa ]  KaI<2  [Va.\  n2  [VA ,]n2 

[Sit.)  NdI\\  [Vco] ni  a(Vc.]n? 

(6.30) 

where  we  have  written  [e~]  =  n,  [h+]  =  p  and  np  =  nf. 

The  relationships  for  the  anti-structures  are: 

[VAs][AsGa]  r, 

[Vfe.]  01 

(6.31) 

.  [VA.\  _  A'ol  I<al 

[I/Ca]  [ASCa]  NeL2 

(6.32) 

,  l$*L]  a  n2 

[S*C?a]  NELin} 

(6.33) 

where  the  extra  equilibrium  constant  is  absorbed  into  a.  The  use  of  the  neutral  EL2 
concentration  is  appropriate  in  this  equation  since  in  n-type  material  the  EL2  remains 
neutral  because  the  Fermi  level  is  much  higher  than  the  EL2  energy  level.  The  con¬ 
centration  of  neutral  EI/2  remains  relatively  constant  and  is  essentially  independent 
of  the  electron  density. 

YYc  know  that  carbon  and  EI/2  are  present  in  significant  amounts  in  GaAs.  The 
relationships  between  ionized  impurities  and  total  impurities  are  given  by: 
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N+=  N * 

l+gexp(£^) 

For  n-typed  semiconductors  at  a  temperature  of  around  1000K,  the  temperature 
where  implanted  dopants  are  usually  annealed,  virtually  all  the  carbon  is  ionized.  On 
the  other  hand,  EL2.  being  a  mid-gap  deep  donor,  less  than  1%  of  it  is  ionized.  As  a 
result,  we  neglect  the  concentration  of  charged  EL2  and  use  the  total  concentration 
of  carbon  in  the  neutrality  equation.  This  condition  breaks  down  however,  when 
the  donor  level  in  GaAs  falls  below  that  of  the  acceptor  level.  In  this  case,  the 
neutrality  condition  must  include  the  ionized  EL2  concentration  and  the  neutral  EL2 
concentration  becomes  highly  influenced  by  the  carrier  density.  Our  ratio  equation 
for  Si  on  the  As  and  Ga  sites  must  also  be  revised  to  include  this  effect.  However,  for 
the  case  of  the  donor  level  being  less  than  that  of  the  acceptor  is  not  of  interest  to  us 
since  this  occures  well  into  the  substrate  not  in  the  active  channel  region. 

If  we  further  assume  that  virtually  all  of  the  Si  dopant  atoms  are  ionized,  then  we 
arrive  at  the  three  equations  which  governs  the  activation  of  the  Si  atoms. 


fo'D  _  <*  ”2 

(S*c«]  Neliti f 

n  +  Nc  +  [SrA.]  =  £  +  [Si&J 

+  [S*Ca]  =  ISiluui 

Letting  z  =  y  =  [Si£J  t  =  [5i]lolai  and  ft  =  a /Nel,2,  resulting  in: 

r? ? 

==>  x  -  y  =  —  —  n  —  Nc  and  x  +  y  =  t 
n 


(6.36) 

(6.37) 

(6.38) 

(6.39) 


1  n2 

^  =  ^i  +  ~-n-Nc) 
1  7 if 

V  =  ^(1  -  ~  +n  +  Nc) 


X 

U 


l  +  t  -  »  -  Ac 
/  -  £  +  n  +  Nc 


=  fin2 


(6.40) 

(6.41) 


(6.42) 
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Multiplying  out  the  above  equations,  we  obtain: 

fin4  +  (l(Nc  +  0*3  +  (1  “  -  <)"?«  -  n •  =  0  (6.43) 

For  each  fitting  constant  fi,  the  above  equation  is  solved  to  obtain  n. 
tuting  n  into  equation  6.36,  we  obtain  the  activation  efficiency  as: 

154.1  -  jgjj  n?  ~  0"’ 

’  ”  [54.1  +  (54.1  »?  + 

where  n,  is  calculated  at  the  annealing  temperature. 

6.2  Se  Activation 

The  first  comprehensive  theory  for  Se  activation  in  GaAs  was  proposed  by  Lidow  et 
al  [35].  They  recognized  the  importance  of  the  complex  formed  by  a  Se  atom  and  a 
gallium  vacancy.  It  is  in  fact  the  mechanism  by  which  Se  diffuses.  It  also  plays  an 
important  role  in  determining  the  amount  of  Se  occupying  As  sites  and  thereby  being 
a  donor. 

The  reaction  process  for  Se  activation  is  described  by  the  following  equations: 


•u  +  Vca  —  (^e^jVco) 

(6.45) 

Sc  A,  =  Se$s  +  e" 

(6.46) 

r  [SeAsVca] 

'  [St-.aV'Ga] 

(6.47) 

1X7  -  [5e.,l 

(6.48) 

[54,1  A-,  (VC,.) 

\SeAaVCa}  KsKi  n 

(6.49) 

where  we  have  used  the  relationship  between  neutral  vacancies.  In  addition,  we  also 


By  substi- 

(6.44) 
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have  the  occurence  of  the  Schottky  anti-site  defect: 


As  A,  **  ASGa  +  2  Va, 


==>  [Vam]  — 


We  finally  arrive  at  the  desired  equation: 


_  K2S/Ksa  1  _  fl'n, 

[SeA,VGa]  I<sl<  i  VNeuti  n 

where  we  have  inserted  an  n,  in  order  to  keep  the  fitting  constant 
neutrality  condition  in  this  case  reads: 


[«•]  +  \CZ\  =  [*+]  +  [Se+At] 


In  additon,  the  total  concentration  of  Se  is  given  by: 


[Self,]  +  [5e4JVc0]  =  [5e]totaj 

=>  (SeJJ  =  n  +  Nc  -  — 

n 

{SeAsVGa}  =  t  -  n  -  Nc  +  — 

n 

where  t  is  the  total  amount  of  Sc  implanted.  As  a  result: 

»  +  Ac  -  =  fTn, 

I  -  n  -  i\(  +  ^  w 

Multiplying  out  the  above  equation,  results  in: 


n3  +  (Nc  fi'n,)n2  -  (n?  +  fi'n.t  -  fi'ntWc)n  -  O'n?  = 


(6.50) 

(6.51) 

(6.52) 
unitless.  The 

(6.53) 

(6.54) 

(6.55) 

(6.56) 

(6.57) 

(6.58) 


Once  again,  the  above  equation  is  solved  to  obtain  n.  Substituting  this  value  of  n 
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Figure  6.1:  Si  Activation  as  a  Function  of  Concentration. 


into  equation  6.51,  we  obtain  the  activation  efficiency  for  Se  in  GaAs: 

_ 

[S<4,]  +  [Se^Vco]  H'n.  +  n 


(6.59) 


6.3  Discussion 

The  activation  of  Si  and  Se  in  GaAs  was  comprehensively  studied  by  J.  K.  Dhiman  and 
K.  L.  Wang  [36].  Their  result  is  implemented  in  the  SUPREM3.5  process  modeling 
program  for  GaAs[37].  By  fitting  Q  and  D'  as  a  constant,  we  may  obtain  the  activation 
efficiency  as  a  function  of  the  total  dopant  concentration.  Figure  6.1  and  6.2  plots 
activat  ion  curves  calculated  using  different  fitting  constants.  The  activation  of  both  Si 
and  Se  is  seen  to  decrease  as  the  dopant  concentration  increases.  Larger  Si  activations 
are  obtained  for  smaller  fitting  constants;  while  larger  Se  activation  is  obtained  for 
larger  fitting  constants.  The  shape  of  our  curves  are  consistent  with  Dhiman ’s  data. 

We  used  Suprem3.5  to  calculate  the  dopant  profiles  for  the  Lot  1  implant  process. 
We  also  fitt  our  activation  mode!  to  calculations  produced  by  Suprem3.5,  whose 
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Figure  6.2:  Se  Activation  as  a  Function  of  Concentration. 


activation  model  is  experimentally  based.  Figure  6.3  and  6.4  show  the  result  of  thes 
calculations.  The  constants  ft  and  ft'  are  selected  so  that  at  the  highest  dopant 
concentrations,  the  activation  efficiency  of  our  model  matches  those  calculated  by 
Suprem3.5.  Figure  6.3  shows  that  the  activation  of  Si  predicted  by  our  model  is  lower 
than  those  that  predicted  by  Suprem3.5.  Figure  6.4  shows  the  activation  predicted  for 
Se  by  our  model  is  higher  than  those  predicted  by  Suprem3.5.  This  result  is  internally 
consistent  with  our  model  and  shows  that  the  site  occupation  of  5*go  >s  proportional 
to  EL2  while  the  Se  occupation  of  the  As  site  is  inversely  proportional  to  the  square 
root  of  EL2.  This  result  is  also  consistent  with  the  reported  depth  distribution  of 
EL2  in  ion  implanted  LEC  GaAs  [38,  39). 

For  these  situations.  theEl/2  concentration  is  typically  higher  at  the  surface,  grad¬ 
ually  falling  to  a  constant  value  well  into  the  substrate.  Since  we  fit  our  ft  and  ft' 
at  the  highest  concentration  of  our  dopant  profile  which  is  always  near  the  surface. 
Deeper  in  the  substrate,  where  the  EL2  level  decreases,  the  fitting  constant  should 
actually  be  larger.  For  larger  ft,  the  Si  activation  is  smaller  which  is  consistent  with 
the  results  presented  in  Figure  6.3.  For  larger  values  of  ft',  the  Se  activation  is  larger 
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Figure  6.3:  Fitting  of  Si  Activation  to  Suprem3.5  Model.  The  value  at  the  highest 
dopant  density  is  fitted  exactly. 


Lot  1  Sc  Activation 


Figure  6.4:  Fitting  of  So  Art i  tion  to  Suprem3.5  Model.  The  value  at  the  highest 
dopant  density  is  fitted  exactly. 
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Fitted  EL2  Concentration 


Figure  6.5:  EL2  Concentration  Resulting  From  the  Fitting  of  Model  Activation  to 
Suprcm3.5  Model. 


which  also  is  consistent  with  the  results  in  Figure  6.4. 

By  varying  the  EL2  concentration,  we  may  fit  the  Suprem3.5  activation  calcula¬ 
tions  exactly.  Figure  6.5  shows  the  result- ng  calculated  EL2  concentrations.  These 
data  are  again  consistent  with  the  experimental  EL2  profiles  obtained  in  [38,  39]. 
Although  the  profile  obtained  as  a  result  of  resulting  from  the  fitting  of  the  Si  activa¬ 
tion  does  not  match  that,  resulting  from  Se  activation,  they  do  bear  a  very  remarkable 
resemblance.  The  final  constant  value  for  the  EL2  concentration  reached  deep  in  the 
wafer  obtained  from  the  Si  and  Se  activation  profiles  also  match  quite  well. 


Chapter  7 

Short  Channel  Effects 

The  short  channel  effect  results  in  a  negative  shift  in  the  threshold  voltage  and  low¬ 
ering  of  the  transconductance,  causing  degradation  of  MESFET  performance.  This 
effect  is  due  fundamentally  to  the  relationship  between  channel  current  and  gate 
length.  It  can  be  seen  even  in  the  simplest  constant  mobility  model  of  the  MESFET 
in  section  2.5.2.  However,  there  are  many  more  reasons  for  the  short  channel  effect. 
Traditionally  this  effect  in  GaAs  MESFETs  is  attributed  to  the  gate  to  source  and 
gate  to  drain  resistances.  The  scaling  of  the  gate  lengths  do  not  proportionally  scale 
the  source  and  drain  resistances.  As  the  gate  lengths  becomes  very  small,  these  re¬ 
sistances  can  have  a  dominant  effect  on  device  performance  [40].  With  the  advent  of 
self  aligned  gate  technology,  these  resistance  are  much  reduced.  However,  the  roles 
of  Carbon  and  EL2  in  short  channel  effect  has  yet  been  fully  explored. 


7.1  The  Role  of  Carbon  and  EL2 

Let’s  first  consider  the  simple  P-N  junction.  The  conduction  and  valence  band  dia¬ 
gram  is  shown  in  Figure  7.1.  The  built  in  voltage  is  given  by: 


v;.  =  e9  -  (vn  +  vp) 


(7.1) 
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EC 


Figure  7.1:  Fermi  Level  with  Respect  to  the  Conduction  and  Valence  Band  for  Simple 
PN  junction. 


Let  the  depletion  region  width  in  the  n-type  side  be  xn  and  p-type  side  be  xv.  Then 
x„  and  xp  and  given  by:  [7] 


2  sVhNa  * 

q(ND  +  NA)ND 


(7.2) 


The  maximum  electric  field  which  occurs  at  the  junction  is  given  by: 


£m 


qNpxn 


qNAxp 

£ 


(7.3) 


This  field  prevents  the  further  drift  of  electrons  from  the  n  region  into  the  p  region. 

Now  suppose  we  add  a  deep  donor  to  the  p  region  whose  concentration  is  above 
that  of  the  p-type  dopant.  The  region  will  now  be  semi-insulating.  The  conduction 
and  valence  band  diagram  for  this  case  is  shown  in  Figure  7.2.  For  this  case,  the 
electrons  will  flow  from  the  n  region  into  the  semi-insulating  regh  this  process, 

some  of  the  ionized  deep  donors  become  neutral  leaving  behind  a  net  ^arge  of  ionized 
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NEUTRAL  N-TYPE 


Figure  7.2:  Fermi  Level  with  Respect  to  the  Conduction  and  Valence  Band  for  Junc¬ 
tion  between  n  type  and  semi-insulating  region. 


acceptors.  In  this  case,  however  the  net  NA  is  much  smaller  than  in  the  p-n  junction 
case.  For  this  situation,  Vj,  is  reduced  bv  approximately  half  the  band  gap.  This 
causes  the  depletion  width  on  the  n  si<  ell  as  the  maximum  field  strength  to  be 

reduced. 

This  configuration  is  identical  to  the  channel  to  substrate  junction.  In  the  presence 
of  EL2,  the  depletion  width  is  reduced  when  compared  to  the  case  of  p-type  substrate. 
This  effectively  increases  the  channel  width.  In  addition,  the  reduced  electric  field 
at  the  junction  makes  it  easier  for  energetic  electrons  to  penetrate  into  the  substrate 
and  form  a  conduction  path.  Both  of  these  factors  contribute  to  an  increase  in  the 
drain  current,  which  in  turn  increases  the  threshold  voltage  of  the  device. 

Figure  7.3  is  a  plot  of  threshold  voltage  as  a  function  of  EL2  concentration.  Here 
the  threshold  voltages  are  calculated  using  the  Medici  simulation  program,  as  are  all 
subsequent  electrical  parameters.  We  see  that  the  magnitude  of  the  threshold  voltage 
is  only  slightly  dependent  on  the  concentration  of  EL2,  which  is  to  be  expected. 
Since  EL2  concentration  above  the  donor  level  is  mostly  neutral,  these  levels  do  not 
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Figure  7.3:  Threshold  voltage  as  a  function  of  EL2  concentration.  The  threshold 
voltages  here  is  derived  from  Medici  simulation. 


contribute  to  current  conduction.  Only  a  small  amount  of  conduction  electrons  arise 
from  the  EL2  donor  increasing  only  slightly  with  increasing  EL2  concentration. 

Figure  7.4  shows  threshold  voltage  calculated  as  a  function  of  carbon  concentra¬ 
tion.  The  concentration  of  carbon  is  seen  to  efTect  the  threshold  voltage  more  than 
EL2  does.  The  concentration  of  carbon  determines  the  metallurgical  junction  of  the 
channel  substrate  junction.  Higher  carbon  concentrations,  the  width  of  the  chan¬ 
nel  is  reduced  which  results  in  a  reduction  of  channel  current  and  threshold  voltage. 
However  as  we  see  in  Figures  4.1  and  4.2,  the  implanted  dopant  profile  is  extremely 
steep  on  the  substrate  side.  So  the  channel  width  is  only  slightly  dependent  upon  the 
carbon  concentration  for  A’e  near  5  x  15  cm-3.  This  result  is  evident  from  inspecting 
figure  7.4. 

Figure  7.5  shows  the  1-V  characteristics  of  a  MESFET  calculated  with  and  without 
EL2.  It  is  evident  that  the  presence  of  EL2  greatly  increases  the  amount  of  current. 
As  we  discussed  earlier,  this  effect  is  due  to  the  reduction  of  the  depletion  width  on  the 
channel  side  when  EL2  is  present.  In  addition,  current  conduction  in  the  substrate 
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Figure  7.4:  Threshold  voltage  as  a  function  of  carbon  concentration.  The  threshold 
voltages  here  is  derived  from  Medici  simulation. 


is  also  increased  with  EL2  present.  This  which  results  from  the  presence  of  EL2 
reducing  the  potential  barrier  at  the  channel  substrate  interface,  results  in  increased 
conduction  current  in  the  substrate.  Figure  7.6  and  7.7,  show  vector  diagrams  of  the 
conduction  current  in  the  channel  as  well  as  substrate,  illustrate  the  above  points. 

7.2  Short  Channel  Effect 

The  short  channel  efTect  for  self  aligned  gate  structures  is  intimately  connected  to  the 
amount  of  carbon  and  EL2  in  the  substrate  region.  As  the  gate  length  is  shortened, 
the  field  in  the  channel  is  greatly  increased  with  the  same  voltages  applied  to  the  gate 
and  drain.  The  increased  electric  field  causes  increased  injection  of  current  into  the 
substrate  region,  resulting  in  non-saturation  of  currents  and  leads  to  a  more  negative 
threshold  voltage. 

Figure  7.8  plots  the  experimentally  measured  threshold  voltage  vs  gate  length 
for  Lot  1  Process  structure.  The  predicted  threshold  voltages  derived  from  Medici 
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Figure  7.5:  I-V  characteristics  of  MESFET  with  and  without  deep  donor  EL2  in  the 
substrate. 


Lot  1  Si  l.lum  with  EL2  -  J.CONDUC 


Distance  (Microns) 


Figure  7.6:  Vector  Diagram  of  Conduction  Current  in  a  MESFET  device.  EL2  is  set 
at  1  x  1016  cm-3  and  Carbon  level  is  at  2  x  1015  cm-3. 
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Lot  1  Si  l.lum  without  EL2  -  J.CCNDUC 


Distance  (Microns) 

Figure  7.7:  Vector  Diagram  of  Conduction  Current  in  a  MESFET  device.  Carbon 
level  is  at  2  x  1015  cm"3.  EL2  is  concentration  is  set  to  0. 

simulation  is  also  shown.  The  gates  of  these  devices  are  oriented  in  the  [010]  direction, 
which  is  45  degress  from  the  major  flat.  As  discussed  in  the  next  section,  these 
directions  are  non-  piezoelectric;  therefore,  can  be  used  in  study  of  short  channel 
effects.  The  Medici  output  shows  good  agreement  with  the  measured  data. 

Using  the  Medici  simulator.  Figure  7.9  and  7.10  are  vector  diagrams  of  the  currents 
in  our  MESFET  devices  calculated.  It  is  evident  that  at  shorter  channel  lengths,  the 
field  between  source  and  drain  is  increased  for  the  same  drain  to  source  bias.  This 
increased  field  causes  increased  current  injection  into  the  substrate.  As  a  result,  the 
current  fails  to  saturate  at  larger  Vga  for  short  channel  devices. 

As  discussed  in  section  7.2,  current  injection  into  the  substrate  is  reduced  if  the 
potent  ial  barrier  at  the  channel-substrate  interface  is  increased.  This  can  be  accom¬ 
plished  by  increasing  the  p-type  concentration  in  the  substrate  by  ion  implantation 
of  a  p-type  dopant  deep  into  the  wafer  near  the  channel-substrate  interface. 
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Threshold  Voltage  Vs  Gate  Length 

La  1  Wafer  114  Si  Imptat 


Figure  7.8:  Threshold  voltage  as  a  function  of  gate  length  for  MESFET  devices. 
Medici  simulation  shows  good  agreement  with  experimental  data. 


Lot  1  Si  0.5um  with  EL 2  -  J.CONDUC 


Figure  7.9:  Vector  diagram  of  conduction  current  for  short  gate  MESFET. 
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Figure  7.10:  Vector  diagram  of  conduction  current  for  long  gate  MESFET 


Chapter  8 

Piezoelectric  Effect 


It  has  long  been  known  that  the  threshold  voltage  of  a  GaAs  MESFET  is  dependent 
on  the  crystalline  orientation  of  the  gates.  This  is  a  problem  unique  to  compound 
semiconductors.  Because  of  the  lack  of  centrosymmetry,  when  external  pressure  is 
applied  to  the  lattice,  ionic  charge  Jistortation  occurs.  This  gives  rise  to  a  polarization 
vector  which  effectively  adds  a  charge  component  to  the  channel. 


8.1  Point  Force  Approximation 

When  thin  films  are  deposited  on  the  surface  of  GaAs  during  device  fabrication, 
eg.  Si3N4,  a  stress  a,  results  in  the  GaAs  substrate.  The  stress  is  given  by  a,  — 
Aajdj/d,  [41],  where  aj  is  the  overlayer  dielectric  film  stress,  and  d j  and  d ,  are  the 
thickness  o'"  the  dielectric  and  semiconductor  wafer,  respectively.  Since  ds  is  ordinarily 
much  greater  than  dj,  the  stress  in  the  GaAs  is  much  lower  than  that  in  the  dielectric. 
However,  wb~n  windows  are  opened  up  in  the  dielectric,  such  as  for  deposition  of  the 
Schottky  gate  metal,  the  stress  imparted  locally  to  the  GaAs  substrate  from  the 
edge  of  the  dielectric  film  is  greatly  enhanced,  approaching  the  stress  value  in  the 
dielectric  [41,  42]. 
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Figure  8.1:  Schematic  diagram  of  a  GaAs  FET  used  in  the  calculation  of  the  stress 
distribution  resulting  from  an  opening  in  dielectric  overlayer. 

8.1.1  Stress  Distribution 

Consider  the  case  shown  in  Figure  8.1,  where  a  window  is  opened  in  the  SbN.*  to 
deposit  the  gate  metal.  We  wish  to  calculate  the  stress  in  the  active  channel  resulting 
from  the  tangential  stress  component  applied  at  the  edge  of  the  gate.  Only  the  case 
of  a  wide  gate  FET  is  considered.  In  this  approximation,  the  stress  tensor  does  not 
depend  on  z,  the  FET  width  direction,  by  symmetry.  This  approximation  is  known  as 
the  “Plane  Stress”  case  and  can  be  solved  analytically  [43].  In  the  coordinate  system 
defined  in  Figure  8.1,  only  aXI ,  <7VV,  cr:.  and  axy  =  ayx  are  non-zero. 

The  stress  tensor  for  the  “Plane  Stress”  case  may  be  solved  as  follows.  If  N(t) 
and  T(t)  are  the  normal  and  tangential  stress  applied  to  the  boundary,  where  t  is 
measured  along  the  boundary,  t  hen: 

orr  +  eryy  =  4/?($(z))  (8.1) 


<?yy  ~  <Jrr  +  2ier,y  =  2{z$'(z)  +  ^(z)) 


(8.2) 
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Where  z  is  a  complex  variable  given  by  x+iy  and  $  and  ^  are  given  by: 

1  f00  N  -iT  .. 


1  roc 

9  =  — - — :  /  — - - < 

2VI  7-oo  t  —  2 

T  1  [<*>  N  +  iT.  1  r00  N-iTAji 

y  =  _  /  - +  —  /  - - -tdt 

2zrt  7-oo  t  —  z  2ffz  7-oo  (<  —  z)2 


The  following  symbols  are  defined  for  our  calculation:  Xi  =  x  —  a,  x2  =  x  -f  a, 
j  =  x  -f  iy,  zj  =  Xi  +  iy  and  z2  =  *a  +  *y.  In  our  case,  N(t)  =  0  and  !T(f)  =  cr  at 
x  =  ±a  and  !T(<)  =  0  otherwise. 


1  r00  T 

•(*)  =  7T  7 - 

2ff  7-00  t  —  Z 


a  r-a+s  di  ra+s  dt 

2k  J-a-S  i  —  Z  la-6  t  —  2 


Tlie  integrals  are  evaluated  by  keeping  aS  —  /=constant  and  lettering  6 
that  limit,  we  obtain: 

»(--)  =  -(---) 

K  \2i  Z?J 

4»(r)  =  -2 <J>(z)  -  zV{z) 

and 

_4/„  (x]  x\\ 

_4fa(*\y7  17 y2\ 

.  ,  Afu  ( x i  x2\ 

=  „„  =  ii  fe?  -  2il) 

1  \  ri  r2  / 

where  i/  is  the  Poisson’s  ratio  and  lias  a  value  of  0.23  for  GaAs  [41]  . 


0.  In 


(8.10) 


(8.11) 
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Figure  8.2:  Schematic  diagram  illustrating  the  orientation  of  the  coordinate  system 
used  in  this  calculation  and  that  of  the  MESFETs. 

8.1.2  Effective  Charge  Calculation 

We  shall  now  transform  our- coordinate  system  for  the  specific  calculation  for  (100) 
oriented  wafers.  We  shall  use  a  coordinate  system  with  the  y-axis  oriented  along 


the  gale  length  and  x-axis  pointing  from  the  gate  downward  into  the  substrate.  The 
coordinate  configuration  is  shown  in  Fig.  8.2.  In  this  coordinate  system,  the  stress 

tensor  is  given  by: 

1/a  / 2/iar2  yix2\ 

=  °rr  =  [  A  4 

*  V  rl  r2  ) 

(8.12) 

ifd  fy'i  I/?  ^ 

(8.13) 

•*./>  (y\  y?\ 

°3  =  <7.-.-  =  2  2 

7T  \r\  r\) 

(8.14) 

_  _  _  _  _  _  -ifd  (xy2i  *y2\ 

6  ~  ry  7T  [r*  r\) 

(8.15) 
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The  piezoelectric  tensor  is  a  tensor  of  rank  three  which  relates  the  polarization 
vector  to  the  stress  tensor.  Since  its  second  and  third  indices  are  symmetric,  it  may 
be  reduced  to  a  two  index  form  according  to  the  following  rule:  [44] 

tensor  notation  11  22  33  23,32  31,13  12,21  , 

(8.16) 

matrix  notation  1  2  3  4  5  6 

In  this  new  matrix  notation,  we  define  d,j=dimm  and  djj=2d,mn  where  m^n.  Both 
mm  and  mn  above  transform  to  j  according  to  our  contraction  rule.  Using  the  same 
method,  the  second  ranked  stress  tensor  is  transformed  to  the  one  index  notation 
with  cr,  =  o}k-  With  these  definitions,  the  polarization  vector  may  finally  be  written 

as: 

P,  =  X>.^  (8.17) 

j=i 

Our  wafers  are  all  (100)  oriented  with  their  major  flat  on  the  (Oil)  surface.  The 
devices  are  fabricated  at  0,  +45,  -45  anti  90  degrees  to  the  major  flat,  corresponding 
to  the  [Oil],  [010],  [Oil]  and  [001]  directions  as  shown  in  Figure  8.2.  For  the  x'y'z' 
coordinate  system  shown,  the  only  non-zero  components  of  the  piezoelectric  tensor 
are  di.|=d25=d36=d.  In  computing  the  polarization  vector,  we  set  up  our  coordinate 
system  such  that  the  y-axis  is  parallel  to  the  direction  of  the  gate  length  and  the 
x-axis  points  downward  from  the  gate  into  the  substrate.  Since  this  is  the  coordinate 
system  in  which  our  stress  tensors  are  expressed,  we  must,  therefore  transform  our 
piezoelectric  tensor  into  this  systei  •». 

(1)  [01 1]  oriented  devices  parallel,  to  the  major  flat. 

The  new  coordinate  system  will  be  x=[100],  y=[011]/\/2,  z=[0lI]/\/2. 

So  the  transformation  matrix  is: 

/  1  0  0 

”=  0  -75  75 

V°  ~7i  -ji 
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The  piezoelectric  tensor  transforms  as  a  tensor  of  rank  three  and  is  given  by: 


dijk  —  atlajmakn<?lmn 


(8.19) 


In  this  coordinate  system,  we  obtain: 


d\2  =  d\22  —  — d/2,  d\ 3  =  di33  =  d/2 


d2 6  =  2d2j2  =  — d,  dss  =  2d3i3  —  d 


The  resulting  polarization  vector  is: 


Pi  =  -(<?3  ~  <72) 


=  ^(^-^>-(4-4)) 

»  V  ri  r2  r«  ) 
P,  =  -dae  =  _  fd) 

7T  V  rl  r2  / 


P3  =  dcrs  =  0 


(8.20) 


(8.21) 

(8.22) 


The  polarization  can  be  turned  into  an  effective  polarization  charge  a s  p  = 
—  V  ■  P.  We  finally  obtain: 

Af^d  /2.tj/i  ,  2.T7/2  6xt/|  xy2  \ 

"  =  — — +  ~ "T  +  1,(T'T)J  (8'23) 

(2)  [Oil]  oriented  devices  perpendicular  to  the  major  flat. 

Here  the  new  coordinate  system  is  x=[100],  y=[011]/v/2,  z=[0ll]/>/2-  The 
resulting  transformation  matrix  is: 


1  0  0 

0  75  75 

0  ~7i  75 


(8.24) 
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In  this  coordinate  system,  we  obtain: 

d\2  —  dijj  =  —d/2,  dia  =  d\23  —  —d/2 
d?6  =  2^212  —  d,  d35  =  2d3i3  =  — d 

As  we  can  see,  these  piezoelectric  tensor  components  have  exactly  the  same 
magnitude  as  in  the  previous  case,  but  have  the  opposite  sign.  Therefore,  the 
effective  polarization  charge  will  have  the  opposite  sign. 

(3)  [010]  oriented  devices  at  45  degree  to  the  major  flat. 

The  new  coordinate  system  is  x=[100],  y=[001],  z=[0T0],  and  the  resulting 
transformation  matrix  is: 

1  0  0  > 

0  0  1  (8.25) 

0-10; 

As  a  result,  dn=d25=d3<3  =  -d.  So  only  the  z-component  of  the  polarization 
vector  is  non-zero,  given  by  -dffg.  Since  all  the  stress  tensors  are  independent 
of  z,  the  effective  charge  will  be  zero. 

The  case  where  devices  are  oriented  at  -45  degree  to  the  major  flat  can  be 
similarly  calculated  and  the  effective  polarization  charge  is  zero  as  expected. 

8.2  Gaussian  Force  Approximation 

As  we  see  in  the  point  force  approximat  ion,  the  effective  charge  distribution  is  singular 
at  the  points  where  the  forces  are  exerted.  This  is  to  be  expected  for  a  point  force. 
Because  of  singularities;  however,  the  point  force  approximation  is  not  suitable  for 
device  simulation.  In  addition,  in  real  situation,  any  force  must  be  exerted  over  a 
finite  extent.  Therefore,  we  shall  use  a  Gaussian  approximation  where  the  forces  are 
centered  at  the  edge  of  the  window  opening  and  fall  off  as  a  Gaussian  towards  both 
sides. 
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We  shall  once  again  start  with  equation  8.3  and  8.4. 
*<J)  = 

J_  r°°  T(t){t  -  x  4-  iy)^ 

~  2ic  J- *  ( t  -  x)2  +  y2 


(8.26) 


(8.27) 

(8.28) 


4f  r^dt  ~  -  **'(*) 

-2*(x)  -  c*'(x) 


(8.29) 


*'(*)  = 


-  rj: 


L  r  JniLdl 

!7T  J-oc  (t  —  x)2 

L  r  ^ w*  ~ ~  + ~  J)y) 

br.L<.  ((*-x)2  +  y2)2 


(8.30) 


\  \  r®  r(t)(2(t-j)y) 


(8.31) 

(8.32) 


With  the  above  setup,  we  obtain: 


auy  =  >uhu  (<!>') 


(J  xi  =  4  R<  (Q>)  —  2y/m($') 
aiy  =  cryj.  =  — 2y/^c($,)  -  2/m($) 
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&zz  —  ^(^ii  ±  Oyy) 

where  u  is  0.23  for  GaAs. 

To  calculate  the  effective  charge  concentration,  we  again  re-orient  our  coordinates 
as  shown  in  Figure  S.2.  The  procedure  is  exactly  the  same  as  section  8.1.  For  the 
[Oil]  oriented  FET,  we  obtain: 


p{x,y)  =  M 


6j (t  -  y)3 
((t  -  y)2  +  x’)* 


(v  +  2)x(t  -y)\  . 

at-y)>  +  x>)>) 


where  M  is  a  fitting  constant.  T(t)  here  is  a  general  force  and  will  be  treated  as 
Gaussian  in  subsequent  simulation. 

It’s  worth  noting  that  if  the  point  force  approximation  is  used  in  any  of  the  general 
formula  derived  above,  we  arrive  at  exactly  the  same  result  as  obtained  in  section  8.1. 


8.3  Discussion 

We  write  the  Gaussian  force  as: 

T’=//ex p[-(i±^)2  —  // exp  [— (~~)2]  (8.33) 

where  the  constant  H  can  of  course  be  absorbed  into  the  overall  fitting  constant. 

Figure  S.3,  S.4  and  8.5  shows  the  piezoelectric  charge  distribution  in  the  channel 
of  the  MESFET  obtained  by  numerical  integration.  The  graph  are  symmetric  with 
respect  to  the  center  of  the  gate  as  expected.  The  charges  are  highly  peaked  around 
the  peak  of  the  force. 

Figures  8.6  and  8.7  constructed  from  experimental  data  taken  from  Lot  2  confirm 
our  calculations.  The  threshold  voltages  here  are  wafer  scale  averages.  We  notice  that 
the  Y'f  increased  for  the  ±45  degree  oriented  FETs  have  nearly  identical  threshold 
voltages,  while  the  0  and  90  degree  oriented  FETs  are  split  on  the  opposite  sides  of 
the  4o  degree  FETs.  These  results  confirm  that  ±45  degree  orientation  are  the  non¬ 
piezoelectric  orientations  and  that  the  0  and  90  degree  orientations  exhibit  opposite 
piezoelectric  effects.  The  gate  length  dependence  of  the  45  degree  FETs  result  solely 
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charge 


Figure  8.3:  Charge  distribution  in  the  channel  due  to  piezoelectric  effect.  Gaussian 
force  approximation  for  1.1  /im  device,  a  —  0.01 /un. 


charge 


Figure  8.4:  Charge  distribution  in  the  channel  due  to  piezoelectric  effect.  Gaussian 
force  approximation!  for  1.1/un  device,  a  =  0.1  /im. 


CHAPTER  8.  PIEZOELECTRIC  EFFECT 


83 


0  0.2  0.4  0.6  0.8  1  1.2  1.4  1.6 

depth  (pm) 


charge 


Figure  8.5:  Charge  distribution  in  the  channel  due  to  piezoelectric  effect.  Gaussian 
force  approximation  for  5.0//m  device,  a  =  0.1/rm. 


by  short  channel  effects  while  a  combination  of  short  channel  and  piezoelectric  effects 
are  responsible  for  the  0  and  90  degree  oriented  FET  gate  length  dependence. 
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Threshold  Voltage  Vs  Gate  Length 

Lot 2  Wafer? SiH  Implant 


Figure  8.6:  Piezoelectric  shift  of  the  threshold  voltage  for  Si  channel  implanted  de¬ 
vices. 


Threshold  Voltage  Vs  Gate  Length 

Lot  2  Wafer  242  Se  Implant 


Figure  8.7:  Piezoelectric  shift  of  the  threshold  voltage  for  Se  channel  implanted  de¬ 
vices. 


Chapter  9 
Conclusions 


The  purpose  of  this  study  is  to  understand  the  electrical  properties  of  GaAs  MESFETs 
in  relationship  to  the  physical  structure  of  the  device  and  the  defect  structure  of  the 
substrate.  We  focused  on  three  issues.  Firstly,  the  equilibrium  properties  of  GaAs 
containing  implanted  species  such  as  Si  and  Se.  We  took  into  account  in  this  analysis 
the  As  anti-site  defect  known  as  EL2  and  developed  models  for  the  activation  of  Si  and 
Se  in  GaAs.  Secondly,  we  studied  the  effect  of  short  channels  on  GaAs  MESFETs.  We 
found  again  that  EL2  and  carbon  doping  play  a  crucial  role.  Thirdly,  we  investigated 
the  orientation  dependence  of  the  electric  characteristics  of  GaAs  MESFETs  known 
as  the  piezoelectric  effect.  We  were  able  to  derive  an  equation  for  Jin  effective  charge 
based  on  the  principles  of  elasticity. 

Dopant  activation  of  Si  and  Se  implants  in  GaAs  is  an  extremely  important  issue. 
Understanding  the  activation  processes  well  go  a  long  way  to  explain  reasons  for 
the  threshold  voltage  fluctuations  in  GaAs  MESFETs.  Manufacturing  yields  are  a 
problem  which  plagues  the  GaAs  IC  industry.  Low  yields  greatly  increase  the  cost 
of  integrated  circuits.  By  understanding  and  improving  threshold  voltage  uniformity, 
we  can  not  only  improve  yields,  but  also  increase  the  level  and  scale  of  integration. 

With  a  combination  of  experimental  measurement  and  the  theoretical  models 
developed  here,  we  can  greatly  improve  existing  process  modeling  software  to  derive 
the  amount  of  activation  as  well  as  understand  the  physical  mechanisms  controlling 
activation.  By  understanding  these  mechanisms,  we  improve  our  ability  to  develop 
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processes  that  meet  the  demand  for  fabrication  of  GaAs  ICs  of  increased  sophistication 
and  complexity. 

In  our  study,  we  point  out  the  shortcomings  of  previous  activation  models.  In  our 
derivation,  we  show  that  EL2,  the  deep  donor,  plays  a  crucial  role  in  the  activation 
of  both  Si  as  well  as  Se  in  GaAs.  We  believe  that  this  model  is  consistent  with  all 
reported  threshold  voltage  fluctuation  data  as  well  as  EL2  distribution  mappings. 

Short  channel  effects  are  another  issue  effecting  the  development  processes  for 
large  scale  integration  processes  for  GaAs  ICs.  For  submicron  gate  length  devices, 
the  short  channel  effect  seriously  degrades  the  performance  of  MESFETs  and  must 
be  understood  and  dealt  with  if  high  performance  and  large  scale  integration  are  to 
be  achieved.  We  show  through  detailed  simulation  that  the  ion  implantation  of  a 
p-type  buffer  layer  is  useful  in  reducing  short  channel  effects. 

The  piezoelectric  effect  is  unique  to  compound  semi-conductors.  We  derive  an 
equation  based  on  the  principles  of  elasticity  to  account  for  this  effect.  By  using  a 
Gaussian  force  centered  around  the  edge  of  the  gate  dielectric,  the  piezoelectric  effect 
can  be  effectively  modeled  using  a  two  dimensional  simulator. 

We  believe  further  studies  are  needed  especially  in  the  area  of  Si  and  Se  activation 
in  GaAs.  We  need  to  systematically  understand  the  distribution  of  various  point 
defects  in  GaAs.  Experimental  efforts  are  needed  to  map  the  lateral  as  well  as  depth 
profile  of  these  defects. 


Appendix  A 
Stress  Calculation 


Let’s  follow  the  notation  outlined  in  chapter  8.  We  will  assume  a  point  force  applied 
at  the  edge  of  the  gate  as  shown  in  Figure  8.1. 

1  r~  N  -  iT  . 


<t>  =  —  — — :  f 

2  m  J- 


2m  J-oo  t  -  z 

1  T 

=  —  / - dt 

2n  J-oo  t  —  z 

a  T  r-“+s  dt  ra+s  dt 
2 TT  \  J-a-6  t  —  Z  Ja-S  t  —  Z 


—  7T  Pn(n  —  f>  +  -)  —  in(<7  +  6  +  z)  —  ln(a  +  <5  —  z)  +  ln(a  —  S  —  z)] 


By  neglecting  high  order  of  6.  we  obtain: 


2z  \Z,  ZoJ  7T  \2t  Zo/ 


(A.l) 


~i  =  riexp(z0i)  z2  =  r2exp(:02) 


(A.2) 


<7xx  +  (Jyy  =  1  Re[<l>{ -)]  (A. 3) 

,/PXP(-?,g>)  exp(-ig2)\ 

?r  \  r,  r2  / 
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z$\z)  +  *(*) 


1  r°°  N  4*  iT 

*  =  -£</-«  TTT" -*«-**'«  <A4> 

=  -2  $(z)-z$'(z) 


=  -2*(z)  +  (z  -  z)V(z)  (A.5) 

7T  \  rt  r2  )  x  \  z\) 

2 fd  [exp(-t0,)  exp(-i02) 


•  1 


+»0 


exp(-2i0,)  exp(  — 


-2  ie2)\ 

rl  ) 


] 


(A.6) 


Iin[S*'(r)  +  #(2)) 
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sin  0i 

sin  02  ,  ( 

cos  201 

cos  201  \  1 
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+  y 

- _i  + - 

r, 

r2  V 

ri 

*2  /J 

2/d  I  sin  g)  -4-  sin  0\  cos  2  0\  sin  02  +  sin  02  cos  202^ 


r,  r2 

sin  gj  •  2  cos2  0j  sin  02  •  2  cos2  02 


*  V 

t( 

Vd  (y^l  _  yj2\ 

*  V  r?  7  2  / 


r2 


) 


(A. 7) 


(A.8) 


^y»  -  °rr  =  2Rc[=4>'(c)  +  ¥(*)] 

=  2Re[-24>(5) +  (*-*)*'(*)]  (A.9) 

Oyy  +  Orx  =  4Re[$(:j)]  (A. 10) 

=»  2ayv  =  2 Re[[s  -  *)$'(*)] 
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4/rf  / ii(r?  -  y2)  x2(rj  -  y2)\ 

"  ar  V  r?  r\  ) 

AJi(i jLiT\ 

7T  V  r?  »*2  y 

=  7  (ffn  4"  °yv) 

4/</7  +  x,y2  Xj  4  *2y2\ 

~  *  {  >1  r<  ) 

Aful  /fi  _ 

~  r  U  r?; 

So  in  summary,  the  stress  distribution  in  the  channel  due  to  point  forces 
at  the  edge  of  the  gate  is  given  by: 
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